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SUMMARY 

_ . research project was to test and truth some recently developed 

— - * *• 

(VAB on me —ere «« - 

re"! IZT^^oZ Sour. Mass Spectrometer and the Neutral 
Atmosphere Temperature Experiment on the satellite. 

Tto activities performed under the one year performance period of the grant may be 
summarized as follows: 

m A major survey of the Atmosphere Explorer V dam base was firs. 
£fold in order to identify the orbits for which suitable 
Experiment, Open Source Mass Spectrometer and Neutral 
Temperature Experiment observations existed. 

(2) Satellite versions of the twilight airglow inversion program -Wchtdiowed 
for the viewing geometry of the VAE observations were generated and tested 
using synthetic data. The inversion program was also modtfted to allow for die 
analysis of twilight observations which included contributions from regtons of 
space with local solar zenith angles less than 90 degrees. 

(3) The twilight observations made on selected orbits were " andAe 
atomic oxygen densities and thermospheric temperatures recovered from these 
Versions were compared with the densities and temperatures measured at . the 
satellite The O+pP) ionization frequencies, which are also recovered as part o 
me inversion process, were compared with the frequencies deduced usmg othe 

methods. 

The results of the study show that a, bod, low and high 
atomic oxygen densities and thermospheric remperanrres recovere 
are in reasonably good agreement with the in r.m sa, ellrte dam. The O ( P) 
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frequencies recovered for low leve.s of solar activity are in 
on the solar F 10.7 0 UX - 


A full journal article describing the results of this work is 
abstract is being submitted to the American Geophysical Union 
Fall 1992 AGU Meeting . 


in preparation and an 
for presentation at the 
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1. BACKGROUND TO THE RESEARCH ACTIVITIES 

*- 7" rr— 

rrr. 7 

j variations of the thermospheric atomic oxygen densittes. Unfortu y, 

however many aspects of both the long term and short term variations, sue as ose 
u r by geomagnetic storms, are sdll not fully understood. At pmsen. there are no 
mrnm in mbit Providing atomic oxygen dam and diem is, therefore, a well mcognumd 

7 . — * “ — 

“7,777,— 777- -m. — • — 

airglow emission features may provide a great deal of information about the long mrm 
"fir, term variadons in thermospheric temperatures and thermospheric oxygen atom 

densities. 

One of die most promising ground-based thermospheric monitoring techrdques 

proposed during the last few years is based ^ method , 

airglow emission a, 7320 A (Femelly et al. 1991, McDade al ) 

orimnally pioneered by Meriwether er al. (1978) and Noxon and Nonon < 1979 >. 

use twilight airglow emission rate measurements made at low elevation in e 

I li ng or setting sun to determine the oxygen atom densines and thennosphen 

1™! Unfortunately, as originally formulated, this method suffers from the 
temperatures, uu u rTTV fl ^ rh e 0 + ( 2 P) ionization 

limitation that detailed information about the solar E densities 

frequencies at the time of the observations is required in order to recover * 

and temperatures. However, McDade e, al. (1991) have demonstrated that this l.m.tatio 

may be overcome, and valuable information about die OW 

also be obtained, if the 7320 A twilight observations are made in tw 

directions - one a, low elevation towards die sun and the other a, higher elevanon ideal y 

towards die local mnith. Pmliminary ground-based measurements ustng ^ app 

are now underway, however, them am at pmsent no sa, elUte measurements of the oxyg 

atom densities being made and it will not be possible, therefore, to threedy ve y 

Z atomic oxygen densities and thermospheric temperatures recovered us ng this 

technique Folnately, during the Atmosphere Explorer H mission stmultaneou 

' ,< Of the 7320 A twilight airglow emission ( Visible Airglow Experiment) 
measurements of the 7izu a twii g e . Hpnsities 

wem made together with measurements of die thermospheric oxygen atom densines 
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(QP en source Mass Spec 

Atmosphere Temperature m ode they should closely 

(VAE) observations were made mam used t0 , es t the 

resemble ground-based observations and may, ’ results of a 

<W> 7320A twilight inversion procedure, using m c 

study carried out to assess the perfomtance of the twthght P 

Atmosphere Explorer E data base. 
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2. THE TWILIGHT AIRGLOW INVERSION ALGORITHM 

The twilight inversion procedures assessed in this work are discussed in detail by 
McDade cal. (1991) and are only briefly described here. The inversion algortthm ts 
based upon die relatively complete understanding of the OW photochem.stry that has 
emerged from the AE-C and AE-D missions (Rusch et al. 1976, 1977). 

The Cri(2P) ion, responsible for the airglow emission at 7620 and 7330 A, is 
primarily produced under twilight conditions as a result of direct P hotolon “ a “ on 
excitation of atomic oxygen by solar EUV photons with wavelengths less than - 6 , 

0 + hv (X<666 A) — >0 + ( 2 P) + e 


( 1 ) 


The ion may also be produced as a result of photoelectron impact lomration « c “ uo ” 
atomic oxygen, but this source is thought to make only a -10% contriburion under most 
twilight conditions {Torr e, al. 1990). The CripP) ions are lost through the rad.au 

decay process 

( 2 ) 


0 +( 2 p) -» 0 + ( 2 D, 4S ) + hv (7320-30A, 2470A) 


and are primarily quenched at thermospheric altitudes by atomic oxygen and molecular 
nitrogen. 


0+(2P) + o o + + o 

0+(2p) + N 2 0 + + N 2 or 0 + N 2 + 


(3) 

(4) 


Because of this relatively simply photochemistry, the twilight 7320-30 A volume 
emission rate at any point in space defined by the altitude, z. and local solar zenith ang e, 

P, may be expressed as 

V(z,P) = {y x A x P(z,p)} / { A + k 0 [O] z + k N2 [N 2 ] z } 

where k, and k N , are die rate coefficients for quenching of OW by atomic oxygen and 
molecular rutrogen; A is the inverse 

ratio for emission of the ( 2 D 5/2 «- 2 P 3 n.ia) ana u 3/z vivi v _., 9m 

and 7330 A ; and P(z,P) is the altitude and solar zenith angle dependent local ( 

volume production rate. 
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The 0+(2P) volume production rate due to ionization excitation by EUV photons in a 
narrow wavelength interval centered on the wavelength X, is given by 


Qi(z.P) = [0] z x fj x I* x exp[ -Tj(z,P)] 


(6) 


where T, is the optical depth for the radiation at wavelength X, and f, is the fractional 

contribution made by radiation in to interval to the - W) " 

„ .. „„ 0 nn C ai depth. The total 0*PP) volume producuon rate, P(z,|3>, may be 
oblerby summing the Q,(z,P) over aU wavelength intervals that contribute 

sienificandy to I*. This is achieved by binning the solar EUV spectrum into the 
wavelength intervals described by Ton « al, (1979). The contribution that photons in a 
given wavelength interval make to the total ionization frequency is calculated using only 
the spectral shape of a reference solar EUV spectrum. 

As most of the twilight 7320-30 A emission originates from altitudes above -250 km 
( Torr et al., 1990; Fennelly et al., 1991) the atomic oxygen densities may be 
approximated by a single exponential profile. The oxygen atom density at any altitude 
Z then be expressed in terms of two parameters - the atomic oxygen scale height Ho, 

which is determined by the thermospheric temperature, and 1 the • “ ojgen atom 
number density, [O h50 , at an arbitrary reference alutude of 250 km. Similarly, me 

molecular nitrogen density at each altitude may be expressed in terms of 1 he density at a 
reference altitude of 250 km and the nitrogen scale height Hni - (l6/28) x H °- 

In die thermosphere most of the attenuation of the EUV flux is due to absorption by 
atomic oxygen and molecular nitrogen and for exponential O and N 2 profiles, the optic 
depth at wavelength Xj may be obtained from the expression 

Tj(z,P) = {[0] a x H 0 x <Ji° x Ch(P, Ho)) + t[N 2 ] a x H Nj x x Ch(P, H m )) (7) 

where o 0 and are the total O and N 2 cross sections at X,; Ch(P, H) is the grazing 
incidence Chapman function and [0] a and [Nj] a are die O and N 2 number densities at the 

minimum ray height of the grazing solar radiation. 

By integrating the volume emission rates given by equation 5 along the line-of-sight 
corresponding to a particular ground-based or satellite-borne twilight observation, the 
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d 7320 30 A column emission rate may be expressed in terms of the well known 
measured 7320-30 A ^ ^ ^ important , and variable, 

physical quantities appearing m (5) (6) U gntl it is possib le m 

geophysical parameters I*. Ho, [Ota. obMm tions by 

principle to deduce the later four quanttt.es from a senes * 

finding die set of four parameters dtatbest it is re^y only 

however, dte two parameters I* and Pta L consists of a series of 
possible to 7^;” depression angles a. two different elevatmn 

— X le a, low elevation towards the azimuth of the sun and the other towards 

particularly “^^^“H^pendently measured nitrogen densities or 

'^densities from a smndard atmospheric model are substituted for [Nitah 

The probiem of finding the values for die parameters I*. Ho. and " 

reproduce a given set of twilight “^^“^^r^nt-^ion method 

rr: - - - * cr. 

■ , for the set of parameters that optimizes the agreement net 

efficiently searches for th p . . . . f the Y 2 merit function, 

model and a set of observations through minimization of X 
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3. THE VAE OBSERVATIONS AND THE TEST DATA REQUIREMENTS 
3.1 The Basic Data Requirements 

In order to make a meaningful assessment of the inversion procedures described in 
Section 2 it was important to use satellite observations that were obtained with vtewtng 
geometries that were as similar as possible to those that would be used to make the 
ground-based twilight measurements. Ideally, the ground-based observations would 
consist of a series of twilight brightness measurements made over a time interval during 
which the solar depression angle at the ground varied between 5 and 20 degrees 
(McDade el al. ,1991). Approximately one half of these observations would be made at 
an elevation angle of -20" towards the azimuth of the setting, or rising, sun and the other 
half would be made towards the local zenith. The idealized ground-based observing 
geometry for which similar satellite observations were to be found is illustrated in Figure 

la. 

3.2 The Visible Airglow Experiment 7320-30 A Observations 

The Visible Airglow Experiment (VAE) on the Atmosphere Explorer 'E' satellite (AE- 
E) was a filter wheel airglow photometer designed to measure various thermospheric 
emission features during both daytime and nighttime conditions. The photometer had 
two distinct optical channels, a high sensitivity channel with a large field of view (3° h 
cone angle) and a low sensitivity channel with a narrow field of view (3/4° half cone 
angle). The fields of view of the two channels were oriented at 90 degrees from each 
other. The counts from the narrow channel (channel 1) were integrated over a period of 
32 msec and those from the wide channel (channel 2) were integrated over a 125 msec 
interval. The instrument operated in a number of different modes depending on the filter 
wheel position which could be held fixed to continuously monitor a particular airglow 
feature or stepped through a number of different interference filters at various stepping 
rates. Continuous observations of the airglow emission near 7320 A were made m two of 
the eight possible VAE fixed wheel modes. In one of these modes, known as mode 
•73F6' the 7320 A observations were made with the narrow channel of the instrument 
(channel 1); in another fixed wheel mode, mode '55F7, the wide channel (channel 2) was 
used to make the 7320 A observations. 

During normal satellite operation the VAE instrument was oriented so that the narrow 
channel pointed aft of the spacecraft and the wide channel pointed towards the earth. 
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'cartwheel spin modes. In th : s j s referred to as the 

was anti-parallel to ^ saK , litc spi „ angular momentum vector 

'normal' spin mode. In the cartw referred to as the 

was par,,, to the “f ""^1, scanned through ,1 

'inverted spin mode. In ei er p w h e n the satellite passed through 

zenith angles within the orbit, plane. Consequen y, insmi ment made 

die terminator in the spinning mode with an acuve A , d . based 

twilight observation of the 0*(®-*P) anglow in a manner s.m. 
twilight monitoring station as shown in Figure lb. 


3.3 Primary Data Selection Criteria 

Given die nature of die basic data moments mid 
the AE-E VAE instrument the first task to be performed was y 

orbits which satisfied the following primary criteria: 


(1) The satellite had to 
circular, orbit. 


be operating in a spinning 


mode in a low altitude, near 


(2) The VAE instrument had to be operation, as the satellite passed ttaough the 
Zv or dawn terminator and the instrument had to be observmg die 7320 
airglow in either Channel 1 or Channel 2 . 


(3) The Open Source Mass Spectrometer and the 
Temperature Experiment both had to be operating 
atmospheric density and temperature data. 


Neutral Atmosphere 
and providing good 
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4. IDENTIFICATION OF ORBITS SUITABLE FOR TESTING THE 

inversion algorithm 

The fust step in the search for suitable orbits for testing me inversion algorithm^ 

r n le ^ E orbits during which the VAE instrument was obsetvtng the 320 A 
to find the At orbits arc listed in the 

airglow and the satellite was operating in the spin mode. 

Ki 7f Annendix 1 along with some pertinent information identifying the VAE channe 
table of Appendix g f { Qf 2; ^ spin mode 0 f the satellite - either normal 

r "“TmToloff umes of me instrument and me satellite alUtude and local times 
or inverted, me on/off u ^ identifl ed mese orbits me next step 

associated " was operaIi „ g when me satellite passed through 

T— and to confirm mat density and temperature dam from the Open Son rce 
me terminator ^ Weu(ral Almosphe re Temperame Expenmen, (NATE) 

were avtdlable a. mat ttme. This was achieved using me AE dambase — -vadaMe 

Explorer £ ' “ * “ « Ablaut (UA, fries for all of me 

examine e The UA data files contain a summary of the data from all the 

orbits listed in Table Al. meuAun ; nt prval The orogram NEWLIST 

AE instruments wim a dam point for every 15 ““"^"^^p.otted on a 

channels were pointing towards the zenith. By — ^ ^ betwee n 

— -s rrsrrsr-— — - 

an indication of whether me satellite was passing from me days.de 
sunset) or from me nightside to me dayside (sunrise). 

I'cfnri in Table 1 were identified the actual VAE 
cinrp the twilight passes listed in laoie i wc 

Once the tw g F . that the instrument was operating 

° bSe T v °l h t! LtrlTL observadons were of sufficient quality to satisfy our 
normally and to ascertain brightness 

..ric Since the UA flies only contained summary data of the VAh ongn 
needs. . vaf data files had to be inspected- 

measurements in me zenith, the full time resoluuon VAE dam files 
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This was performed using rhe Fortran program 'VAEREAD' which unpacks rhe VAE 
data fries and extracts the VAE brightness measurements and the ancillary tnstrumen 
“ “ A listing of the version of VAEREAD used in dtrs project, 
TWIVAEREAD, is attached as Appendix 2. The program reads the channel 1 or ctonel 
2 photometer counts, converts the observed counts into Rayleighs and makes 
and galactic background corrections where possible. The program creates two output 
text (ASCII) files. One of the files lists sequentially the following quantities. 


(1) The universal time in milliseconds 

(2) The observed 7320 A brightness in captured Rayleigh untts 

(3) The estimated error in the observed brightness based on the photometer count rate 

(4) The satellite altitude in kilometers 

(5) The zenith angle of the photometer line of sight in degrees 

(6) The solar zenith angle at the satellite in degrees 


The second output file contains the following quantities: 


(1) The universal time in milliseconds 

(2) The x, y and z coordinates of the satellite position in the Geocentric 

Equatorial Inertial (GEI) system (see Russell , ) 


(3) The x, y and z coordinates of the sun in the GEI system 

(4) The x, y and z coordinates of the photometer line-of-sight in the GEI system 


Detailed examination of the high tesolution VAE data for the orb.ts ^ Table 1 
revealed drat in many instances, particularly in the case of the channel 1 dam the 
observations were sJously and urevocably contaminated by stars or ft* 
and galactic background corrections could not be performed. Furthermore, althoug 
VAE instrument was equipped with a sophisticated baffle system, close -ammationof 
the data revealed that the twilight observations made at low elevations towar 
were often contaminated by scattered sunlight. Because of these varrous undestmMe 
effects most of the twilight passes listed in Table 1 had to be rejected. Hav.ng tttentif.ed 
the passes which contained potentially usefirl data a final selection errtenon then had 

applied. 
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As described in Section 3 the ground-based thermospheric monitoring technique^ and 

. *i;_ t. t k measurements made in tne 

Ute twiligh, inversion procedure, requtre ' 732 0 A measum^ ^ 

-t ^ ° f f r 

rri,r;;=rr-:r rrjriss:: 

terminator varied considerably and on many of the of rbe 

SC “°T o ri“drcX" l”-f L twilight passes listed in Table 1 

it— r:« — 

observations made on orbit 7012. 
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Table 1. 


AE-E 


spin twilight passes with OSS, NATE and VAE 7320 A observations 


Date 

yyddd 


Orbit Channel Spin 

# # mode 


77003 

77007 

77007 

77007 

77009 

77009 

77009 

77016 

77016 

77016 

77040 

77047 

77047 

77052 

77052 

77063 

77063 

77068 

77077 

77077 

77081 

77081 

77087 

77087 

77100 

77100 

77108 

77108 

77132 

77132 

77132 

77138 

77138 

77142 

77142 

77146 

77155 

77176 

77176 

77176 

77181 

77181 

77364 

77364 


5811 

5874 

5878 

5884 

5912 

5913 
5913 
6019 
6025 
6029 
6415 
6525 
6525 
6597 
6612 

6776 

6777 
6855 
7012 
7012 
7075 
7075 

7168 

7169 

7372 

7373 
7497 
7497 
7895 

7895 

7896 
7980 
7980 
8055 
8055 
8122 
8257 
8600 
8600 
8601 
8672 
8672 

11617 

11618 


2 

1 

2 

1 

2 

2 

2 

2 

1 

2 

2 

2 

2 

1 

2 

1 

1 

2 

1 
1 

2 
2 
1 

1 

2 
2 
1 
1 
1 
1 
1 
2 
2 
1 
1 
2 
1 
1 
1 
1 
2 
2 
2 
2 


normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

normal 

inverted 

normal 

inverted 

normal 

normal 

normal 

normal 

normal 

normal 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

inverted 

normal 

normal 


Start 

time 

End 

time 

Condition 

25725 

27165 

sunrise 

19005 

19830 

sunrise 

40500 

41295 

sunrise 

72585 

73485 

sunrise 

50085 

50910 

sunrise 

53430 

54540 

sunset 

55155 

55455 

sunrise 

18495 

18915 

sunset 

50565 

51165 

sunset 

71955 

72660 

sunset 

68505 

69495 

sunrise 

52785 

53865 

sunset 

54540 

54975 

sunrise 

8910 

9345 

sunrise 

83970 

84315 

sunrise 

19590 

20505 

sunset 

21495 

22410 

sunrise 

12660 

13560 

sunset 

73320 

74235 

sunset 

75225 

76140 

sunrise 

65850 

66810 

sunset 

67725 

68685 

sunrise 

47970 

48870 

sunset 

51405 

52320 

sunrise 

21795 

22455 

sunset 

26790 

27165 

sunset 

3045 

4065 

sunset 

4875 

5895 

sunrise 

60435 

61005 

sunset 

61935 

62895 

sunrise 

65445 

65820 

sunset 

945 

1305 

sunset 

1995 

3075 

sunrise 

59115 

60545 

sunset 

60545 

61980 

sunrise 

76115 

77550 

sunrise 

24075 

24990 

sunset 

50265 

50565 

sunset 

51435 

52440 

sunrise 

54975 

55645 

sunset 

6960 

8375 

sunset 

45045 

45525 

sunrise 

55620 

56820 

sunrise 

59355 

60450 

sunset 
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5. INVERSION CODE MODIFICATIONS FOR INVERTING 

THE SATELLITE DATA 

As already mentioned, the ground-based twilight inversion algorithm is designed to 
deal with observations made in the zenith and a, low elevation in the azimuth of the sun. 
Since the spin plane of the twilight VAE observations was rarely oriented perpendtcular 
,o die plane of the terminator, the inversion algorithm had to be modified to allow for the 
finite angle between die azimuth of the sun and the photometer line-of-s.ght. 

In order to allow for this effect die program which was used to mad the VAE data 
files, TWIVAEREAD, was modified to extract the Geocentric Equatorial Inertial 
coordinates of (i) the satellite position vector, (ii) the sun vector, and (in) the vector o 
the VAE photometer line-of-sight, as well as die observed 7320 A column emission rates 
and photometer zenith angles. These vectors were then used within the inversion 
algorithm to calculate the local solar zenith angle at various intervals along the 
photometer line-of-sight. The section of code dealing with this problem is incorporate 
in the line integral calculation performed by the procedure 'BRIGHT which appears on 
pages 5 and 6 of the listing given in Appendix 3. 

Other modifications had to be made to die inversion code to allow for the fact that 
many of the VAE observations were made at smaller solar depression angles than wou 
be accessible from the ground. In the case of ground-based observations ^spheric 
sprin g of sunlight makes it very difficult to obtain good CrfPD-2P> 7320-30 A 
measurements until the solar depression angle at the observing site is greater than about 7 
degrees for zenith observations and about 12 degrees for low elevation angle sunward 
observations (.Meriwether et at. ,1978; Noxon and Norton. 1979; Fennelly et al. 199 , 
McDade et at.. 1991). However, tropospheric scattering does not interfere with the 
satellite measurements and good twilight data were obtained for satellite solar depression 
angles down to zero degrees. It was considered valuable to include these 'early twilight 
observations but die early observations made a, low photometer elevation angles 
inevitably included contributions from regions of space lying on the days.de of the 
terminator where the local solar zenith angle, P, was less than 90 degrees. Hie lute 
integral calculation section of the original inversion code, procedure BRIGHT, was 
therefore modified to deal with this situation by including brightness contnbunons from 
both sides of the terminator. For contributing elements with local solar zenith angles 
greater than 90 degrees the grazing incidence Chapman function was used to calculate 
the optical depth as explained in Section 2 (equation 7). For elements with loc so 
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zenith angles less than 90 degrees the normal Chapman function was used and the optical 
depth at each wavelength, Xj , was calculated from the expression 

Xi(z,|3) = { [0] z x Ho x Gi° x Ch(P, Ho) } + l [NJz x H Nz x °i Nz x Ch ^’ Hn2 ^ 

where z is the altitude of the contributing element and [0] z and [N 2 1 Z are the inferred 
atomic oxygen and molecular nitrogen densities at that height. The relevant section o 
the inversion code appears on pages 5 and 6 of Appendix 3. 
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6. RESULTS FROM THE INVERSIONS OF SELECTED SATELLITE ORBITS 

In an ideal world it would be more desirable to test the inversion algorithm using the 
observations made with the narrow channel of the VAE instrument (c anne • 
However, the narrow channel was about 50 times less sensitive titan the wide t channe 
and the advantages of using the channel with the smaller field of view (and short r 
integration period) were strongly outweighed by the much h.gher signal to no.se ra 
the channel 2 observations. Nevertheless, it was considered msm.ct.ve 
algorithm using both the wide channel and the narrow channel observations. 

6. 1 Wide Channel Results - Orbit 6855 

For the purposes of this study the best VAE channel 2 twilight observations were the 
sunset observations made on orbit 6855 on day 68 of 1977. When the measurements 
were made the satellite was in a near circular orbit at a latitude of 12 °N, a on 8 ltu c ° 
215 oe and an altitude of 257 km. The solar F,o.t flux value for the day was 80.3 and the 
Ap index was 38. The angle between the orbital plane and the terminator was 
approximately 80 degrees as the satellite crossed the terminator and the angle bem-een 
the azimuth of die photometer scan and the azimuth of the sun was 10.2 degme 
Although it was important to have the latter angle as small as possible (see Section 4) 
Tmean tha, the observations made a, the lowest elevations in the sunward direction 
were prone to contamination by scattered sunlight. To avoid the danger of : using 
observations with solar contamination we only considered the observations which were 
made when the angle between the sun and photometer line-of-sight was greater an 
degrees. As a result, the analysis was restricted to zenith observations and sunward 
observations made at elevation angles greater than -35 degrees. In .gure we s o 
part of the sequence of 0+PD-2P) 7320 A column emission rates measured g 
sunset pass of orbit 6855. The solid symbols show tire 7320 A emission rates measurj 
when die photometer elevation angle a. the midpoint of die sample mtegration penod 
was within ± 5 degmes of the zenith direction; die open symbols show the “ rate 
measured in the sunward direction when die photometer elevation angle was between 35 
and 45 degrees. Both sets of observations are plotted against the solar zenith ang e a 
satellite which was increasing with time. The satellite was spinning a. approximately 
revolutions per minute and die measurements shown in Figure 2 span a total observing 

period of about 5 minutes. 
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A comparison of the data shown in Figures 2 and 3 reveals the stmtlanty between the 
VAE twilight measurements and the simulated ground-based observations discusse y 
McDade a al., (1991). It should be noted, however, that relative to the zenith 
observations the sunwatd observations of Figure 2 are weaker than those of Figure 3. 
This is primarily due to die fact that the sunward VAE observations were made a. an 
elevation angle of 40±5 degrees and the sunward simulations were calculated for an 
elevation of only 20 degrees. It is also evident that the uncerrainties in the VAE 
measurements are considerably larger than those of die ground-based simulations m spite 
of the fact drat the satellite observations were made in the absence of any tropospheric 
background scattering. However, it has to be recognized that each of die satelhm 
observations was obtained with a sample integration period of only 0.125 seconds 
whereas the ground-based simulations were calculated for a one minute integration 
period using an instrument with a larger throughput. I. should also be pointed out that 
not all of die scatter evident in the VAE observations is to be associated with noise 
because much of it is due to variations in the photometer elevation angle within die 
plotted ±5 degree elevation angle bands. This occurs because the satellite was spinning 
at approximately 24 degrees per second and the photometer counts were sampled every 
0125 seconds, therefore, there were usually three observations made on each spin 
between the elevation angles of 35 and 45 degrees and three observations made within ±5 
degrees of the zenith. Since the observed twilight brightness depends on both the 
photometer elevation angle and the solar depression angle, the ±5 degree spread in 
elevation angles contributes to the scatter which should be greater in the case of the 
sunward viewing observations. 

When the 7320 A emission rates of Figure 2 were corrected for the On®- J P) 7320 
A and 7330 A doublet filter capture functions and processed using the twilight inversion 
program (see Appendix 3) the algorithm returned the following set of fitting parameters t 

(i) an atomic oxygen scale height of Ho = 62 ± 3.4 km 

(ii) an oxygen atom density at 250 km of (Ol a0 - 1.4 ± 0.3 x 10» cm-3 

and (in) an unattenuated 0 + PP) ionization frequency of I* = 6.2 ± 0.6 x HP 8 sec-' 


The fitting parameters were not found to depend significantly upon the first guess 
values used to initiate the inversion procedure although the number of iterations 
required to reach convergence did of course depend upon the initial guess, me fits to 
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ft. zenith and sunwa* observations obtained using these best ^parameters are shown 
in Figures 4 and 5. In the case of the sunward observations. Figure 5, the fit s 
toothed' rather than smooth because of the ±5 degree spread in the elevation angles 
discussed above. Since the twilight brightness should no. ^so snong.ybe^een^S 
and -5 degrees of the zenith this effect does not show up m the t 
observations. 

Clearly the parameters recovered from the fit to die orbit 6855 observations do 
reproduce fairly well die input data. More significandy, however, they also reproduce 
die orbit 6855 observations that were no. used in die inversion, r.e. the ^eredfi 8 
parameters also reproduce the VAE observations that were no. used to obuun dm 
parameters. This is illustiated in Figures 6, 7, 8 and 9 which show ““ “ * 
sunward observations that were made a. elevation angles in the ranges 45 to 55 degrees, 

55 to 65 degrees, 65 to 75 degrees and 75 to 85 degrees. 

The neutial atmospheric temperatures measured a. 257 km by dm Neutral 
Atmosphere Temperature Experiment (NATE) as die satellite passed through die sunse 
terminator on orbit 6855 varied between about 870 K and 950 K. Since die recovered 
atomic oxygen scale height of 62 ± 3.4 km is equivalent to an exosphe^emfmrafrrm » 
970 ± 50 K we see that the thermospheric temperature inferred from the twi g 
observations is in very good agreement with the NATE temperamre —mem. 
Similarly, the average atomic oxygen density measured by the Open Source 
Spectrometer (OSS) during the twilight pass was 9 x 10 s cm’ at 257 km w ic comp 
very favorable with the inferred density of 1.2 ± 0.3 x I0» cur 3 based on the mcovered 
density at 250 km and a scale heigh, of 62 km. In Figure 10 we show how dm atomic 
oxygen density profile constructed from the recovered fitting parameters Ho and [O te „ 

compares with die profile based on dm NATE temperamre and the OSS density 
measurements. Figure 10 also shows the atomic oxygen denstty profile given by th 
MSIS-86 model (Hedin, 1987) for the sunset conditions on orbit 6855. 

The value for the unattenuated 0*GP) ionization frequency obtained from the 
inversion of the orbit 6855 channel 2 observations wiU be discussed m Section 7. 
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6.2 Narrow Channel Results - Orbit 7012 

The best channel 1 twilight 7320 k observations were those obtained on the sunset 
nass of orbit 7012 on day number 77 of 1977. When these observations were made the 
satellite was in a circular orbit at a latitude of 16 °N, a longitude of 322 “E and an 
altitude of 251 km. The solar F 10 ., flux value for the day was 75, the Ap index was 
and the angle between the azimuth of the photometer scan plane and the azimuth of die 
sun was 14.5 degrees. As already mentioned the signal to noise ratios of the channel 1 
observations were very much lower than those of the channel 2 observations. Th.s is 
clearly fflustiated in Figure 11 which shows the orbit 7012 channel 1 observauons made 
within ±5 degrees of the zenith direction. Because of the low signal to noise ratios in e 
data it was not possible to simply invert the channel 1 observations made within ±5 
degrees of the zenith and between 35 and 45 degrees elevation towards the sun, t.e. the 
very large uncertainties associated with the recovered fitting parameters rendered the 
inversion meaningless. However, somewhat more meaningful results were obtaine 
when all of the observations acquired between elevations of 35 degrees and the zenith 
were considered. The entire set of observations obtained between sateBtte solar zenith 
angles of 90 and 105 degrees in the zenith and sunward in the elevation angle an 
45 , 45-55 , 55-65 , 65-75 and 75-85 degrees are shown in Figure 12. 

When this entire set of observations was inverted the inversion algorithm returned the 
following set of best fit parameters. 


(i) an atomic oxygen scale height of Ho - 76 ± 14 km 

(ii) an oxygen atom density at 250 km of [Ol^o = 70 ± ®- 4 x 108 cm 

and (Hi) an unattenuated CCPPj ionization frequency of I* = 9.6 ± 3.5 x 10-8 sec' 

The fit to the entire set of orbit 7012 observations obtained using these parameters is 
shown by the solid line through the smoothed data points in Figure 12. 

The average neutral atmospheric temperature measured at 251 km by the Neutral 
Atmosphere Temperature Experiment ( NATE) as the satellite passed through the sunset 
terminator on orbit 7012 was 880 K and the average atomic oxygen density measured y 
the Open Source Mass Spectrometer (OSS) during the same period was 1.0 x 10 cm . 
Clearly, the recovered oxygen density of 7.0 ± 0.4 x 1(>8 ctn-3 at 250 km agrees wtth the 
OSS measured density within the uncertainty limits. The recovered atomic oxygen scale 
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height of 76 ± 14 km is equivalent to an exospheric temperature of 1 180 ± 220 K which 
is somewhat hotter than the NATE measured temperature of 880 K. 


6.3 Results at Higher Levels of Solar Activity - Orbit 2 4564 

In order to rigorously assess the twilight inversion procedures it was considered 
highly desirable to test the inversion algorithm using VAE observations made at both ow 
and high levels of solar activtty. Unfortunately, most of the twilight passes wh,ch 
satisfied the primary selection criteria outlined in Section 3, and which contained data of 
a sufftciendy high quality, wete made at low levels of activity. However, there were a 
small number of passes made in 1980 close to solar maximum for which oxygen atom 
data were not available but which did involve channel 2 observations made when the 
satellite was spinning. One of these high activity orbits which contained apparen y goo 
data was orbit number 24564 on day 1 00 of 1980. The twilight obstrvations made on the 
sunset pass of orbit 24564 towards the zenith and towards the sun at an elevauon of 0 
degrees are shown in Figure 13. When these observations were made the AE-E sate 
wm, in a circular orbit at a latitude of 6.6 °N, a longitude of 80.4 °E and an alutude of 
419 km. The solar F,„., flux value for the day was 244, the Ap index was 20, the ang e 
between the azimuth of the photometer scan plane and the azimuth of the sun was 3.5 
degrees and the satellite was operating in the inverted spin mode. 

When the orbit 24564 column emission rates shown in Figure 13 were inverted the 
inversion procedure returned the following set of best fit parameters: 


(i) an atomic oxygen scale height of Ho - 94 ± 8 km 

(ii) an oxygen atom density at 250 km of [O] 25 o = 2-5 ± 0.3 x 10 cm 

and (iii) an unattenuated 0*( 2 P) ionization frequency of I* = 9.3 ± 1.0 x 10' 8 sec' 1 

For the inversion of this data the shape of the EUV flux spectrum used in the tdgonUtm 
(see Appendix page 14, was based on the 79050 spectrum teported by Ton et al (1979) 
however, the results obtained from inversions using the shape of the Stan ar 
spectrum of Hinteregger (1977) were not signiftcandy different. The fits to the zem 
and sunward observations obtained using the best fit parameters listed above are shown 

in Figures 14 and 15. 
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ZZZTZ O K and «. — s 

1470 ± 120 K inferred from the recovered atom,c orb ® bot *«. recovered 

, ^cc Awopn atom data was not available for tms orDii uui u 
Unfortunate * 2 50 ^ compares very weU with the density of 2.2 x 

conditions by the MS1S-86 model (Hedin, 1987, Howeve, « w 
10 cm , . thp recovered 0 + ( 2 P) ionization frequency of I* - 9.3 ± 

will discuss m the next se » might be expected for conditions close to solar 

1.0 x 10-8 se c-i is substantially lower than might be expectea 


maximum. 
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7. DISCUSSION AND CONCLUSIONS 

me results described in the previous section clearly demonstrate that the Aversion 
procedures for recovering thermospheric temperatures and atomic oxygen densities from 
directional ground-based measurements of die OW) 7320 A riv, light -glow 
emisrion petfood well when tested with die proxy satellite data. The atomic oxygen 
densities recovered from the inversions are in reasonable 

densities measured by the Open Source Mass pectrome er heights are 

thermospheric temperatures infetred from die recovered atomic oxygen scale heighri am 
also in reasonably good agreement with the measurements made on dm sarellde by * 
Hemal Amosphere Temperature Experiment. Furthermore, the temperatures are also 
good agreement wid, die temperatures that have been deduced using o*er rechnt^ 
This is illustrated in Figure 16 which shows how the temperature, , recove :re 
6855, 7012 and 24564 compare with the temperatures deduce y eea 
from an analysis of the late twilight 7320 A zenith intensities measured on these an 

other AE-E orbits. 

•me unattenuated 0+(*P) ionization frequencies recovered for orbits 6855 and 7012 
are in good agreement with previous evaluations but the ionization frequency recovered 
for orbit 24564 is somewhat smaller than might be expected for conditions close to solar 
“Hum. me solar cycle dependence of the 0* P P) ionization frequencies has been 
studied by Ton et al. (1979) who used the solar EUV flux measurements on e 

Atmosphere Explorer satellites W »") » " H 7 1980 
frequencies on five selected days during the 1974 to 1979 period. Abreu et al. (1980) 

have also investigated die solar cycle dependence of the W ""77 

and used dayglow 7320 A measurements made with the VAEmstrumen o 

Leniencies during die increasing phase of solar cycle 21. The ionization frequencies 
HlTfrom di! work of Torr et til. (1979) and Abreu e, tri. (1980) are shown tn 

Figures .7 and ,8 where they are comped 

VAE twilight observations on orbits 6855, 7012 an 

recovered from orbits 6855 and 7012 are consistent with what should be expecte at ow 
levels of activity but the orbit 24564 frequency is not in keeping with e tten s 
previous evaluations. Because of the lack of appropriately conditioned high ac vity 
Light observations it is difficult to determine whether or not the seeming y ow 
ZLa frequency for orbit 24564 is indicative of a problem wid, a*. — 
algorithm. It is important to note, however, that die atomic oxygen scale height and 
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densities recover^ from the inversion of the orbit 24564 data are h, good ^men, 
with the temperatures measured on the satellite and the densities pre c 
86 model. We should also point out that if the orbit 24565 observations are inverted with 
the 0<-(2P) ionization frequency constrained to a value that is in keeping wi e tren 

in Figure 18. ihen die recovered atomic oxygen scale heigh, and densities are no 

longer in good agreement with die measured and modelled ^ die 

For example, if the ionization frequency is constrained to 1.3 x 
inversion algorithm returns an oxygen atom density at 250 km of 2.5 x 1 Cm “ 
thermospheric temperature of only 1320 K. I, is possible, however, tha, the 

low value for the orbit 24564 ionization frequency is simply a reflecaono () 

Zrint’of die orbit 24564 frequency from die trend tine in Figum^lS . not 
inconsistent with the scatter in the measured frequencies at low io.i 
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no. 1. (a) Sketch illustrating 
etal. (1991). 



nG . (h) Same as Fig. la hut musing how «h, observations were made 
with the Visible Airglow Experiment on the AE-E satel . 
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FIG. 2. The zenith and sunward twilight 7320 A emission rates 

VAE channel 2 on the AE-E sunset pass of o , < >855^ The cmission 

r ang.es ranging from 35 to 45 

degrees. 


30 




Final Report NASA Grant NAG 5-1502 



no 3 The simulated ground-based twilight 0*PP) 7320 A observations d^sed by 
VcD™e e,T(1991). The sunward column emission rates were calculated for an 

elevation of 20 degrees towards the azimuth of the setting sun. 
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FIG. 4. The channel 2 VAE 7320 A column emission rates measured within ±5 degrees 
of the zenith on AE-E orbit 6855 (data points) and the fit obtained using the 
parameters discussed in the text (solid line). 
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FIG. 5. The channel 2 VAE 7320 A column emission rates measured in the sunward 
direction at elevation angles between 35 and 45 degrees on orbit 6855 (data points) 
and the fit obtained using the parameters discussed in the text (solid line). 
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FIG. 6. The channel 2 VAE 7320 A column emission rates measured in *e sunward 
direction at elevation angles between 75 and 85 degrees on orbit 6855 (da a _p ) 

and the fit obtained using the parameters discussed in the text (solid line). . 


observations were not used to find the fitting parameters. 
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SOLAR ZENITH ANGLE (deg) 


FIG 7 Hie channel 2 VAE 7320 A column emission rates measured in the sunward 
n °dJecn?n afeTevahon ang.es between 65 and 75 degrees °n orbit 6855 (data potnt^ 
and the fit obtained using the parameters discussed in the text (solid line). N.B. these 
observations were not used to find the fitting parameters. 
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mr s The channel 2 VAE 7320 A column emission rates measured in the sunward 
nG d^cto afe—V between 55 and 65 degrees — » 
and the fit obtained using the parameters discussed in the text (solid line). N.B. th 
observations were not used to find the fitting parameters. 
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FIG. 9. The channel 2 VAE 7320 A column emission rates measured in the sunward 
direction at elevation angles between 45 and 55 degrees on orbit 6855 ( ata points) 
and the fit obtained using the parameters discussed in the text (solid line). N.B. these 
observations were not used to find the fitting parameters. 
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FIG. 10. The atomic oxygen density profile (solid line with error bars) recons^cted 
using the H 0 and [O ] 250 parameters obtained from the fit to the orbit 6855 tw g 
observations. The profile derived from the temperature and atomic oxygen densities 
measured by the NATE and OSS instruments is shown by the open square and dashed 
line. The atomic oxygen densities from the MSIS-86 model are represented by the 

dot-dashed curve. 
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PXG. 11. The channel 1 VAE 7320 A column emission rates measured within ±5 degrees 
of the zenith on AE-E orbit 7012 (data points). The emission rates were obtained by 
averaging the channel 1 counts over four 32 msec integration periods. The plotted 
data have been smoothed (in both the ordinate and abscissa) using a three point 
running average to illustrate the underlying trend. 
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FIG. 12. The orbit 7012 channel 1 VAE 7320 A column emission rates measured withm 
±5 degrees of the zenith and sunward in the elevation angle bands 35-45 , 45-55 , 55- 
65 65-75 and 75-85 degrees. Each sequence of 50 points shows how the emission 
rates varied between the solar zenith angles 90 and 105 degrees. The plotte ta 
have been smoothed as in Fig. 11 to illustrate the underlying trend. The solid line 
shows the unsmoothed fit obtained using the parameters discussed in the text. N.B. 
for the inversion described in the text raw unsmoothed data were used. 
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trrr 14 The channel 2 VAE 7320 A column emission rates measured within ±5 degrees 
of the ^l on ^forbU 24564 (data poims) and the ft obtained ustng the 

parameters discussed in the text (solid line). 
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FIG 16 The temperatures recovered in this work for AE-E orbits 6855, 7012 and 24564 
(squares with error bars) compared with those deduced by Tee and Abreu , (1 982) 
from their analysis of zenith 7320 A twilight observations (solid tnangles) 
temperatures measured by the Neutral Atmosphere Temperatu ire E x P en ™ 
shown by the open triangles highlighted with central dots on orbits 6855, 7012 and 


24564. 
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FTG 17 The 0 + (2P) ionization frequencies recovered in this work for AE-E orbits 
fiS55 7012 and 24564 (squares with error bars) compared with the frequencte 
deducecMty Mreuetal. 1980) from dayglow 7320 A observadons (open mangles) 
and”*/ frequencies calculated by Ton e, al. (1979) from EUV flux measurements 


(solid triangles). 
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FTr. 18 The 0 + ( 2 P) ionization frequencies recovered in this work for AE-E orbits 
6855 7012 and 24564 (squares with error bars) compared with the frequencies 
deduced b ^ IXeu et al 1980) from dayglow 7320 k observations (open tnangles 
amHhe^ frequencies calculated by Torr a/. (1979) from EUV flux measurements 


(solid triangles). 
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Table A1 

AE-E orbits with spinning satellite and VAE observations at 7320 A 


Date Orbit* Ch Spin Time (sec ITT) 

off 


yyddd # 

# 

type 

75352 353 

2 

N 

75354 372 

1 

N 

75354 378 

2 

N 

75356 395 

1 

N 

75356 397 

2 

N 

75356 402 

1 

N 

75356 406 

1 

N 

75358 420 

1 

N 

75358 422 

2 

N 

75358 427 

1 

N 

75358 431 

1 

N 

75360 447 

2 

N 

75360 452 

1 

N 

75360 456 

1 

N 

75362 470 

1 

N 

75364 496 

1 

N 

76022 798 

1 

N 

76024 818 

1 

N 

76032 925 

1 

N 

76036 973 

1 

N 

76038 1001 

2 

N 

76040 1026 

2 

N 

76042 1048 

2 

N 

76042 1054 

2 

N 

76048 1125 

2 

N 

76052 1174 

2 

N 

76054 1202 

2 

N 

76093 1721 

2 

N 

76096 1764 

2 

N 

76099 1801 

2 

N 

76099 1810 

2 

N 

76108 1931 

2 

I 

76111 1967 

2 

I 

76114 2008 

2 

I 

76117 2049 

2 

I 

76120 2091 

2 

I 

76123 2133 

2 

I 

76126 2174 

2 

N 

76126 2183 

2 

N 

76132 2269 

2 

I 

76135 2302 

2 

I 


on 


57008 

16080 

57712 

2816 

16672 

51304 

78992 

3112 

16944 

51512 

79216 

17080 

51640 

79656 

3167 

9735 

84205 

45492 

78820 

56588 

72171 

67331 

41835 

82043 

37843 

18387 

31466 

39296 

50136 

25008 

82240 

69063 

36143 

32511 

28815 

29974 

30958 

24230 

79582 

84238 

27158 


59367 

18447 

60071 

5183 

19039 

53663 

81359 

5479 

19311 

53871 

81495 

19439 

53999 

81294 

5534 

12102 

84836 

47851 

81099 

58955 

74530 

69698 

44202 

84410 

40210 

20746 

33825 

41303 

52135 

27015 

84239 

71062 

38150 

34510 

30814 

31981 

32957 

26237 

81581 

85405 

29157 


Altitude (km) 

on mid/ap off 


Local Solar Time (hr) 


952.1 

1042.1 

1053.8 

1023.8 

1022.6 

1027.3 

1042.9 

1051.8 

1066.9 

1097.1 
1059 

1032.5 

1019.6 

587.3 

972.3 

894.4 

386.5 
956 

686.5 

843.2 

769.2 

693.8 

875.5 

851.4 

926.9 

891.9 

830.8 

532.2 

380.2 

543.1 

519.9 

667.7 

317.7 
738 

517.7 

720.9 

477.8 

727.6 

744.5 

380.8 

1146.5 


140.9 
141.2 
140.9 
140.8 
141 
140.8 

140.8 
145.1 
145.1 

145.1 
145 
145 

144.9 

144.7 
145 
145.5 
145 
144.4 

144.1 

144.1 
144.4 

145.1 

140.8 

141.1 

143.3 
142.7 

141.3 

141.9 
141 
139.6 

139.4 

137.1 

137.9 
139.6 

141.4 

141.9 

140.5 

138.2 

137.6 
152.1 
151.8 


1098.9 

1001.1 

987.2 

1029.1 

1029.8 
1023.3 

1006.2 
997 

981.6 

932.9 

898.4 

1009.7 
1003 

641.8 

1045.4 
1136 

1001.6 

985.8 

1150.6 

1054.2 

1106.8 

1195.6 

974.5 

1012.9 

908.6 

913.1 

962.1 

776.7 

1871.3 

723.4 

731.6 

530.5 

927.7 

438.3 

369.1 

433.4 

644.4 

399.8 

373.9 

236.9 

166.6 


on 

mid < 

off 

2.9 

7.88 

13.56 

2.85 

8.21 

13.53 

2.89 

8.34 

13.56 

3.19 

8.6 

13.91 

3.22 

8.64 

13.94 

3.26 

8.78 

13.98 

3.27 

8.86 

13.99 

3.44 

9.06 

14.13 

3.42 

9.1 

14.11 

3.4 

9.21 

14.02 

3.58 

9.28 

13.97 

3.91 

9.48 

14.55 

4.03 

9.6 

14.6 

5.76 

9.64 

13.4 

4.49 

9.81 

15.03 

5.23 

10.17 

15.78 

17.6 

19 

20.35 

10.01 

15.09 

20.66 

12.48 

17.02 

22.8 

12.7 

17.78 

23.29 

13.53 

18.21 

0.01 

14.29 

18.53 

0.83 

13.93 

18.85 

0.5 

14.12 

18.9 

0.75 

14.87 

20.21 

1.61 

15.71 

21.16 

2.36 

16.43 

21.7 

3.01 

2.85 

6.9 

11.95 

10.66 

14 

18.04 

4.15 

8.17 

13.46 

4.4 

8.33 

13.68 

5.66 

10.77 

14.83 

8.63 

11.39 

17.63 

6.83 

12.07 

15.86 

8.67 

12.79 

16.29 

8.26 

13.75 

17.62 

10.37 

14.69 

19.57 

9.7 

15.47 

18.88 

9.79 

15.57 

18.86 

13.76 

17.03 

19.18 

10.11 

17.74 

18.66 


"when an observing sequence spans m«Tta one orbit only fe fust orbi. number is listed 
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APPENDIX 1 


Table Al. continued 


Date Orbit 

yyddd # 


Ch 

# 


Spin 

type 


Time (sec UT) 
off 


on 


76135 

76138 

76141 

76144 

76153 

76153 

76156 

76156 

76159 

76159 

76162 

76162 

76165 

76168 

76168 

76171 

76171 

76177 

76177 

76180 

76180 

76186 

76186 

76189 

76189 

76192 

76192 

76210 

76210 

76213 

76213 

76216 

76216 

76219 

76219 

76222 

76222 

76225 

76225 

76228 

76228 

76231 

76231 

76234 


2311 
2344 
2387 
2426 
2557 
2563 
2594 
2600 
2638 
2644 
2680 
2689 
2721 
2768 
2774 
2809 
2815 
2896 
2902 
2937 
2943 
3020 
3026 
3062 
3068 
3104 
3110 
3358 
3364 
3400 
3406 
3442 
3448 
3485 
3491 
3529 
3535 
3571 
3577 
3613 
3619 
3656 
3662 
3699 


2 

2 

1 

1 

1 

2 

1 

2 

1 

2 

1 

2 

1 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 


82309 

26853 

32325 

12365 

39036 

75812 

6532 

43276 

18500 

55372 

10716 

65916 

4684 

27892 

64819 

21076 

57835 

44043 

80947 

36699 

73539 

28450 

65178 

27898 

64570 

26818 

63658 

27049 

63825 

25625 

62377 

23288 

60001 

27296 

63952 

35936 

72520 

26320 

62872 

21464 

57720 

22632 

58888 

22280 


84316 

28852 

34324 

14364 

41035 

77811 

8235 

44979 

20507 

57371 

12715 

67923 

6691 

29899 

66714 

23075 

59842 

46050 

82946 

38698 

75546 

30449 

67185 

29905 

66577 

28817 

65665 

29048 

65832 

27624 

64384 

24711 

62007 

29295 

65959 

37895 

74519 

28319 

64879 

23471 

59727 

24639 

60887 

24279 


Altitude (km) Local Solar Time (hr) 

on mid/ap off on mid off 

991.5 

151 

758.7 

11.04 

11.6 

12.2 

654.2 

151.7 

450.6 

13.28 

18.65 

22.57 

395.7 

151.5 

725.4 

15.9 

19.49 

0.89 

553.2 

152 

543.5 

15.61 

20.17 

0.62 

458.3 

156.6 

639.5 

18.57 

22.77 

3.84 

460.3 

156.9 

636.9 

18.67 

22.89 

3.92 

923.1 

155.9 

175.3 

16.67 

23.5 

0.49 

1042.2 

155.5 

155.7 

16.16 

23.67 

23.75 

474.3 

153.4 

644.8 

20.21 

0.31 

5.28 

467.9 

153 

639.8 

20.38 

0.38 

5.38 

1107.1 

152.1 

180.9 

17.44 

0.9 

2.08 

1259.7 

151.9 

152.3 

16.74 

1.06 

1.2 

532.8 

151.7 

583.3 

21.37 

1.62 

6.66 

738.4 

151.6 

387.7 

20.95 

2.62 

6.25 

672.4 

151.6 

369.7 

21.39 

2.75 

6.23 

538.7 

151.6 

570.3 

22.81 

3.52 

8.06 

667.1 

151.4 

457.3 

22.16 

3.63 

7.48 

472.7 

153 

653.7 

1.12 

5.01 

10.22 

484.3 

153.3 

629 

1.16 

5.09 

10.21 

657.8 

155.1 

453.8 

0.74 

5.7 

9.83 

678 

155.2 

444.6 

0.73 

5.82 

9.9 

798.7 

157.1 

345.8 

1.4 

7.51 

10.63 

936.4 

156.5 

262.1 

0.8 

7.68 

9.97 

542.4 

154.9 

562.8 

3.66 

8.39 

12.79 

690.5 

154.4 

433.3 

2.87 

8.49 

12.08 

698.8 

153.3 

433.2 

3.62 

9.06 

12.7 

751.5 

152.9 

388.9 

3.44 

9.16 

12.49 

638.4 

153.4 

464 

8.82 

13.79 

17.83 

657.1 

153.8 

446.7 

8.83 

13.9 

17.83 

442.9 

155.6 

653.4 

10.8 

14.55 

20.02 

445.7 

156 

660 

10.88 

14.65 

20.18 

632.4 

156.9 

195.1 

10.33 

15.44 

16.95 

625.1 

157.2 

461.7 

10.47 

15.62 

19.76 

343.7 

145.8 

758.8 

13.04 

16.45 

22.19 

316.7 

145.7 

795 

13.42 

16.61 

22.5 

557.7 

143.5 

467.4 

12.39 

17.34 

21.31 

527.7 

143.2 

524.2 

12.73 

17.44 

21.8 

405.6 

141.3 

644.6 

14.35 

17.99 

23.38 

338.2 

141.3 

736.7 

15 

18.09 

0.1 

658.9 

140.6 

389 

13.46 

18.74 

22.59 

698.7 

140.8 

346.8 

13.33 

18.85 

22.36 

417.8 

140.3 

604.4 

15.67 

19.64 

1.06 

393.6 

140.6 

630.1 

15.95 

19.76 

1.33 

610.8 

140.4 

389.8 

15.12 

20.68 

0.43 
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APPENDIX 1 


Table Al. continued 

Date Orbit Ch Spin Time (sec UT) Altitude (km) 
yyddd # # type on off on mid/ap off 


76234 

3705 

2 

I 

58767 

76237 

3739 

1 

I 

3823 

76243 

3832 

1 

N 

41103 

76243 

3838 

2 

N 

76839 

76246 

3875 

1 

N 

38311 

76249 

3926 

1 

N 

79423 

76255 

4006 

1 

N 

26038 

76255 

4012 

2 

N 

61150 

76264 

4138 

1 

I 

17510 

76264 

4144 

2 

I 

52246 

76271 

4248 

1 

I 

42701 

76274 

4290 

1 

I 

27157 

76274 

4299 

2 

I 

77893 

76277 

4338 

1 

I 

32757 

76277 

4347 

2 

I 

83101 

76283 

4433 

1 

N 

39501 

76283 

4439 

2 

N 

72788 

76286 

4478 

1 

N 

29604 

76286 

4486 

2 

N 

73452 

76292 

4573 

1 

N 

32460 

76292 

4579 

2 

N 

65068 

76295 

4625 

1 

N 

55756 

76298 

4675 

1 

N 

67691 

76301 

4721 

1 

N 

57939 

76304 

4773 

1 

N 

78403 

76319 

5008 

1 

N 

37514 

76319 

5014 

2 

N 

69466 

76322 

5057 

1 

N 

41202 

76322 

5063 

2 

N 

73154 

76340 

5340 

1 

N 

6849 

76346 

5440 

2 

N 

24041 

77003 

5810 

2 

N 

22542 

77007 

5874 

1 

N 

19054 

77007 

5878 

2 

N 

40502 

77007 

5884 

1 

N 

72630 

77009 

5912 

2 

N 

50094 

77016 

6019 

2 

N 

18533 

77016 

6025 

1 

N 

50614 

77016 

6029 

2 

N 

72006 

77034 

6307 

2 

N 

7516 

77034 

6313 

1 

N 

39636 

77034 

6317 

2 

N 

61044 

77036 

6352 

1 

N 

76996 

77040 

6415 

2 

N 

68508 


60390 

378.9 

140.6 

400.8 

5830 

372.8 

139.7 

627.8 

43102 

268.3 

141.4 

707.6 

78846 

245.1 

142 

750.1 

38622 

244.3 

142 

375.1 

81430 

386.4 

142.7 

528.4 

28045 

442.2 

137.8 

438.1 

63157 

488.2 

137.3 

390.7 

19517 

433.4 

135.3 

363.4 

54253 

389.9 

134.7 

394.5 

44700 

447.4 

136 

273.7 

29164 

393.1 

142.3 

297.4 

79900 

428.3 

141.9 

262.8 

34508 

414.6 

141.7 

192.3 

85108 

481.4 

140.6 

197.1 

41508 

349.8 

139.8 

251.6 

74723 

351.2 

139.2 

228.7 

31603 

275.7 

137.9 

281.2 

75451 

342.2 

137.5 

219.7 

34467 

233.2 

142.8 

266.7 

67075 

250.1 

141.8 

246.9 

57763 

273.4 

148.4 

207.7 

69690 

246.1 

161.5 

229.1 

59090 

168.3 

161.4 

211.3 

80402 

258.7 

169.8 

204.1 

39521 

241.8 

184.1 

190.3 

71473 

246 

184.1 

187.2 

43201 

243.3 

238.8 

247.5 

75161 

247.2 

239.1 

244.2 

12152 

250 

247.2 

250.2 

29344 

244.4 

241.6 

244.6 

27901 

249.6 

253.5 

249.6 

19917 

247.8 

246.2 

245.8 

41341 

247.2 

245.7 

245.5 

73493 

247.2 

245.4 

244.8 

55461 

244.5 

247.3 

244.4 

19396 

258.6 

256.9 

255.9 

51477 

259.4 

257.3 

255.8 

72869 

259.8 

257.7 

256.1 

8379 

263.7 

262.1 

260.4 

40499 

263.8 

262.4 

260.6 

61907 

263.4 

262.1 

260.4 

82355 

259 

257.9 

259.2 

72347 

259.5 

254 

255.5 


Local Solar Time (hr) 
on mid off 


16.87 

20.77 

0.62 

17.66 

21.5 

2.78 

20.5 

23.13 

5.67 

20.87 

23.32 

6.09 

2.79 

3.5 

4.19 

21.12 

1.36 

6.41 

22.38 

2.91 

7.48 

22.12 

3.02 

7.21 

0.68 

5.83 

10.05 

1.19 

5.97 

10.47 

2.61 

8.01 

11.69 

3.95 

9.02 

13.24 

3.69 

9.25 

12.97 

4.42 

10.3 

12.65 

3.76 

10.55 

12.99 

7.02 

12.31 

16.09 

7.1 

12.43 

15.81 

8.9 

13.21 

18.12 

7.9 

13.45 

17.1 

11.34 

15.72 

20.45 

10.97 

15.83 

20.16 

11.29 

16.7 

20.27 

13.16 

17.83 

22.35 

17.71 

19.17 

23.16 

14.62 

20.63 

23.82 

18.65 

1.84 

3.73 

18.3 

1.95 

3.31 

23.18 

2.32 

8.25 

22.11 

2.54 

7.1 

0.38 

9.02 

0.14 

20.81 

6.77 

20.6 

9.36 

21.3 

9.3 

3.32 

5.09 

6.96 

3.26 

4.97 

6.78 

3 

4.77 

6.65 

3.21 

15.23 

3.23 

20.58 

22.38 

0.22 

19.81 

21.64 

23.47 

19.3 

21.14 

22.97 

7.91 

9.79 

11.58 

7.15 

9.06 

10.93 

6.62 

8.55 

10.44 

8.86 

20.78 

8.7 

1.62 

10.04 

18.78 
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Table Al. continued 

Date Orbit Ch Spin Time (sec UT) Altitude (km) 

yyddd # # type on off on mid/ap off 


77043 

6455 

1 

77043 

6459 

2 

77047 

6524 

2 

77052 

6598 

1 

77052 

6602 

2 

77052 

6612 

2 

77063 

6776 

1 

77068 

6855 

2 

77077 

7011 

1 

77081 

7074 

2 

77087 

7168 

1 

77108 

7297 

1 

77100 

7372 

2 

77108 

7497 

1 

77132 

7895 

1 

77138 

7980 

2 

77142 

8055 

1 

77146 

8122 

2 

77155 

8257 

1 

77169 

8492 

2 

77176 

8600 

1 

77181 

8672 

2 

77187 

8770 

1 

77194 

8882 

2 

77214 

9211 

1 

77218 

9279 

2 

77224 

9366 

1 

77228 

9436 

2 

77275 

10191 

1 

77275 

10191 

1 

77279 

10249 

2 

77285 

10345 

1 

77314 

10816 

1 

77318 

10880 

2 

77324 

10979 

1 

77328 

11041 

2 

77350 

11396 

1 

77354 

11456 

2 

77360 

11550 

1 

77364 

11617 

2 

78005 

11705 

1 

78054 

12026 

2 

78028 

12085 

1 

78032 

12144 

2 


N 

20716 

21435 

N 

41916 

42779 

N 

49619 

54978 

N 

8483 

9346 

N 

29907 

30770 

N 

83467 

84330 

I 

19554 

24921 

I 

11962 

17329 

N 

73225 

76944 

N 

63705 

69073 

N 

48009 

53376 

I 

2224 

5767 

I 

21824 

27183 

I 

2224 

7591 

I 

60470 

65837 

I 

990 

6349 

I 

59013 

64381 

I 

73821 

79180 

I 

22204 

27571 

I 

74108 

79475 

I 

50307 

55674 

I 

7003 

12370 

I 

17282 

22498 

I 

16242 

19793 

I 

59721 

65088 

I 

79641 

80768 

I 

28704 

34071 

I 

61616 

66975 

N 

65941 

68260 

N 

68269 

68908 

N 

31181 

36452 

N 

31332 

32299 

N 

60643 

66002 

N 

58186 

62881 

N 

75834 

76769 

N 

64298 

68153 

N 

73536 

74127 

N 

48472 

53839 

N 

43752 

45543 

N 

55663 

61030 

N 

11190 

13109 

N 

31466 

33825 

N 

78389 

83748 

N 

54260 

56899 


252 

252.5 

252.9 

251.5 

252.2 

252.9 

247.8 

250.6 

247.7 

244.7 

243.1 

242.5 

244 

242.5 

242.2 

243.7 

242.1 

241.5 

261.7 

261.8 

261.6 

258.3 

259 

258.3 

251.5 

253.9 

250.4 

245.4 

249.6 

245.3 

240 

241.6 

240 

258.8 

258.8 

261 

264.8 

263.2 

264.9 

258.8 

258.4 

258.7 

277.4 

273.3 

277.4 

270.2 

271 

270.2 

269 

267.6 

269.1 

267.9 

264.7 

267.9 

261.5 

261.8 

261.5 

252.8 

251.1 

252.9 

281.6 

282.1 

281.8 

279.3 

280.9 

279.4 

278.3 

276.4 

279.1 

275.8 

274.6 

276.9 

266.7 

263.2 

266.8 

264.6 

263.6 

263.2 

288.6 

282.7 

288.7 

284.8 

281 

284.9 

251.4 

252.2 

253.3 

253.3 

253.5 

254.1 

278.8 

254.8 

280.6 

279.6 

278.2 

278.3 

269.6 

268.1 

270.3 

281.3 

278.6 

279.1 

274.7 

273.6 

273.6 

272.7 

271.3 

278.3 

281 

280.4 

280.3 

279 

280.8 

279.1 

279.9 

280.1 

275.7 

273.9 

273.7 

274.1 

270.3 

268 

267.5 

830.8 

141.3 

962.1 

320.6 

324.8 

321 

320.3 

319.3 

322.4 


Local Solar Time (hr) 
on mid off 


9.81 

11.48 

13.18 

8.5 

10.4 

12.45 

3.73 

15.74 

3.74 

0.55 

2.4 

4.32 

0.44 

2.29 

4.21 

0.24 

2.03 

3.96 

17.54 

5.46 

17.43 

14.65 

2.63 

14.61 

17.22 

1.27 

9.65 

7.95 

19.9 

7.85 

2.52 

14.49 

2.51 

14.3 

22.24 

5.95 

19.3 

7.21 

19.17 

14.3 

2.29 

14.26 

19.16 

7.03 

18.97 

20.62 

8.51 

20.39 

17.23 

5.2 

17.14 

14.29 

2.2 

14.15 

9.71 

21.62 

9.59 

0.66 

12.61 

0.59 

20.34 

8.2 

20.12 

17.95 

5.82 

17.74 

14.28 

1.85 

13.45 

9.67 

17.58 

1.33 

21.47 

9.38 

21.34 

18.56 

20.91 

23.45 

11.35 

23.2 

11.1 

11.05 

22.91 

10.83 

8.34 

13.4 

18.74 

18.81 

20.15 

21.48 

3.63 

15.42 

3.14 

2.29 

4.3 

6.39 

7.97 

19.9 

7.82 

20.76 

7.22 

17.77 

2.77 

4.73 

6.78 

0.27 

8.9 

17.33 

10.89 

12.08 

13.37 

18.34 

6.19 

18.12 

20.32 

0.44 

4.41 

1.82 

13.71 

1.64 

22.69 

2.55 

6.82 

16.43 

21.7 

3.01 

8.53 

20.35 

8.16 

6.61 

12.41 

18.25 
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Table Al. continued 


Date Orbit 

yyddd # 


Ch 

# 


Spin 

type 


79096 
79100 
79102 
79104 
79106 
79108 
79110 
79112 
79114 
79116 
79132 
79134 
79136 
79138 
79142 
79146 
79148 
79150 
79152 
79154 
79156 
79158 
79160 
79162 
79164 
79176 
79178 
79178 
79181 
79184 
79187 
79190 
79193 
79196 
79199 
79202 
79205 
79220 
79223 
79226 
79232 
79235 
79238 
79253 
79256 
79259 


18868 1 
18929 1 

18960 1 

18992 1 

19023 1 

19054 1 

19085 1 

19115 1 

19146 1 

19177 1 

19424 1 

19455 1 

19486 1 

19517 1 

19578 1 

19640 1 

19670 1 

19701 1 

19731 1 

19762 1 

19793 1 

19824 1 

19855 1 

19886 1 

19917 1 

20103 1 

20133 1 

20133 1 

20180 1 
20222 1 
20267 1 

20318 1 

20364 1 

20410 1 

20457 1 

20503 1 

20549 1 

20781 1 

20827 1 

20871 1 

20966 1 

21012 1 
21059 1 

21290 1 

21335 1 

21381 1 


Time (sec UT) 

off 


on 


on 


Altitude (km) 

mid/ap off 


Local Solar Time (hr) 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

N 

N 

N 


74599 
71815 
73247 
74535 
75959 
77391 
78478 
74214 
75550 
76878 
76013 
77261 
78501 
79741 
76620 
79180 
74860 
76164 
71860 
73124 
74371 
75595 
76819 
78019 
79211 
80578 
76106 
78986 
80673 
57154 
50513 
77601 
76633 
75633 
80168 
79024 
77872 
77711 
76511 
64055 
78639 
77510 
81950 
81405 
74565 
78933 


78094 
75350 
76598 
76886 
79702 
79542 
82109 
77837 
79173 
80493 
79636 
80884 
82124 
83364 
80307 
82923 
78603 
79907 
75483 
76755 
77994 
79218 
80442 
81642 
82834 
84201 
78977 
79729 
84297 
60785 
54136 
81224 
80256 
79256 
83791 
82647 
81503 
81334 
80134 
68038 
82261 
81133 
85573 
85028 
78188 
82555 


453.6 

452.6 
452.3 

452.1 

452.6 
453.8 

453.7 

454.7 

455.8 
456.6 

452.2 
452.1 

451.8 

452.3 

452.1 

451.4 

451.2 

451.4 

452.3 

453.3 

454.3 

454.9 
454.9 

454.5 

453.4 
447.1 
447 
453.7 
447.3 
447.1 

447.1 

447.3 

448.1 

450.1 

452.2 

453.4 

453.5 

447.2 
446.9 

446.3 

445.3 

445.4 
446.3 

447.6 

443.7 

443.5 


457.7 

459.8 

459.9 
457.8 

458.4 
458 

456.4 

454.7 

453.4 

452.2 

449.5 

449.7 

454.7 

450.7 
453 

455.2 

455.5 
455 

454.2 
453 
451.9 

451.1 

450.6 

450.2 

450.2 

448.8 

447.8 

455.2 

448.9 
450 

451.5 

451.9 

451.1 

448.8 

446.9 

445.9 

445.5 
446.8 

447.7 

449.8 

452.3 

451.8 

450.2 

448.6 
447 

442.4 


457.6 
454.9 
454.9 

458.7 
454.1 
456 

454.6 

454.7 

454.4 

454.1 

457.1 

456.8 

456.1 
455 
452 

449.5 

449.2 
449.1 
449 
449 

448.8 

448.5 

448.4 

448.6 

449.1 

455.2 

453.7 

455.4 

454.5 

453.2 

451.1 

449.5 

449.2 

449.2 

448.5 

447.8 

447.5 

450.8 
450 
447 

445.4 

445.3 

445.5 

449.8 

451.5 

452.8 


on 

mid 

off 

4.62 

11.95 

19.64 

4.94 

12.29 

19.95 

5.19 

12.15 

19.29 

4.81 

9.64 

14.85 

4.93 

13.05 

20.86 

5.21 

9.85 

14.51 

4.04 

11.86 

19.38 

4.06 

12 

19.69 

4.04 

11.92 

19.73 

4.09 

11.81 

19.71 

3.66 

11.62 

19.2 

3.47 

11.43 

19.1 

13.68 

16.24 

18.95 

3.24 

11.02 

18.88 

3.21 

10.9 

18.98 

3.48 

11.31 

19.33 

3.66 

11.59 

19.49 

3.85 

11.93 

19.75 

4.07 

11.99 

19.56 

4.13 

12.08 

19.81 

4.25 

12.09 

19.87 

4.24 

12.05 

19.94 

4.35 

11.97 

19.93 

4.35 

11.88 

19.84 

4.29 

11.85 

19.75 

3.25 

11.22 

18.83 

3.05 

9.44 

15.37 

15.37 

17.01 

18.7 

3.08 

10.86 

18.72 

3.04 

10.73 

18.68 

3.23 

10.8 

18.72 

3.67 

11.24 

18.95 

4 

11.86 

19.42 

4.31 

12.21 

19.85 

4.37 

12.25 

20.08 

4.27 

11.94 

19.89 

4.16 

11.76 

19.72 

3.53 

11.49 

19.2 

3.38 

11.23 

19.09 

3.17 

11.57 

20.27 

3.8 

11.41 

19.15 

4.03 

11.76 

19.34 

4.1 

11.98 

19.57 

16.31 

17.78 

19.3 

12.2 

15.79 

19.27 

3.91 

11.78 

19.34 
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Appendix 2 

The following is a listing of the Fortran program used to extract 7320 A emission 
rates and related orbital quantities from the AE-E Visible Airglow Experiment data files* 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Program TWIVAEREAD 

a version of the AE data base program VAEREAD . FOR 
created for the project 

-An Assessment of Twilight Airglow Inversion Procedures 
Using Atmosphere Explorer Observations- 

by I.C.McDade 

under NASA Grant NAG 5-1502 

This program extracts selected data and vectors from VAE data files 
The output files are used for graphical examinationofthe data 
for input to the twilight inversion program TWIFITTER . PAS 

Logical unit assignments, all done internally: 

LOGICAL* 1 ID 
CHARACTER *80 FSTR, F2STR, VSTR 
INTEGER *4 MODEX(7,3), LAM {7, 3) 

REAL *4 OADATA (22) ,OADATL<22) ,RAY1(4) ,RAY2,STV(2) ,DC(2) , 

>sun (3 ) , vl ( 3 } , v 2 (3) , vll (3) , v21 (3) , r (3) , sgeo(3) , sgei (3) , 

>sth, sph , Tmat (3,3) , gst 

C sun is the sun position vector in GEI coordinates 

C V 1 YJ 2 are the channel 1 & 2 line-of-sight vectors in GEI coords. 

C r is the satellite spin axis vector in GEI coordinates coords 

c sgeo & sgei are the satellite position vectors in GEO & GEI coords 
C sth & sph are the satellite sherical coords in degrees 
C Tmat is the 3x3 GEO to GEI transformation matrix 

COMMON /CVEC/ VI , V2 , V1L, V2L, R, sun, gst 

COMMON /CRAYL/ ICH, NSAT , IGZC, STV, DC , RAY1,RAY2 .ICE.iSKIP 
COMMON /CVAERD/ ITIME, THET, ICH11 , ICH12 , ISQ11 , ISQ12 , ITIMEL, THETL, 
>TFW IE ICH11L, ICH12L, ISQ11L, ISQ12L, IATN1 , 

>ICH2 , ISQ2 , IATN2 , TBAF1 , TBAF2 , TBAF11 , TBAF22 , TAEL, TBEL , TPMT1 , 

>TPMT2 , TFW, 

>IEND4, IX, ISVIF,OADATA,OADATL , 

DATA MODEX/ ' 55F6 ' , ' 73FD 1 , 1 52F7 1 , ' 63F4 1 , ' 33F5 ' , 'CLF2 , 42F3 , 

> ' 73F4 ' 1 52FD 1 , ' 42F5 ' , 1 48F6 ' , ' 55F7 ' , 1 CLF4 ' , 1 63F5 ' , 

> 1 73F6 ’ 1 52FD ' , '28F5' , '65F6' , '55F7' , ' CLF2 ' , ' 63F5 1 / 

DATA LAM/5577, 7320, 5200, 6300, 3371, 0, 4278, 

>7320, 5200, 4278, 4861, 5577, 0, 6300, 

>7320, 5200, 2800, 6563, 5577, 0, 6300/ 


C Open input and output files 


5 

10 


TYPE 5 

FORMAT ( ' Enter name of 
ACCEPT 10, VSTR 
FORMAT ( A80 ) 

OPEN (UNIT=4, FILE=VSTR 


input VAE data file:') 
TYPE= ' OLD ' , READONLY , 
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+ FORM= 1 UNFORMATTED 1 ) 

PRINT*, 'Enter NREADMAX ?' 

READ*, NREADMAX n 

PRINT*, 'Do you want G&Z cor (1-Y, 0 N) . 

READ*, IGZC 

PRINT*, 'Enter channel # (1 or z) 

READ* , ICH 

prS*; 'Intl/S? w ch2 zenith angle correction (if any) ?■ 
READ*, zacor 

11 PRINT*, 'MIN channel 1 zenith angle . 

READ*, ZMIN 

PRINT*, 'MAX channel 1 zenith angle . 

READ* , ZMAX 

15 FORMAT^ { ' Enter name of output file containing observations 
S P ^ONIT=2 o! R NAME=FSTR , TYPE= ' NEW ' ) 

16 FORMAT 6 (' Enter name of output file containing the vectors :') 
ACCEPT 10, F2STR 

OPEN (UNIT=2 1 , NAME=F2STR, TYPE='NEW') 


C Assign initial values: 


NS=0 

NERR=0 

NSQ=0 

NATN=0 

NAVG=0 

STO=0 . 

ISVIF=1 

ISKIP=0 


C Get start, stop times, and averaging parameter: 


35 

40 


45 


TORMAtV Enter start, stop times in seconds (defaults, 86400) :' ) 
ACCEPT *, JSTART , JSTOP 
ISTART = JSTART*1000 
IF (JSTOP .EQ. 0) JSTOP = 86400 
ISTOP = JSTOP*1000 

FORMAT 5 (' Enter Filter Wheel Position (10=all) : ' ) 

ACCEPT *, JFW 

C Position Vae Data File to start time, obtain header information and 
C check mode 

CALL INITV ( 4 ! ISTART, IDATE, IORB, ID, ISPINF, INF , MODE1 , ITON , 

> ITOFF , IDMF , MODE2 , ITON2 , ITOFF2 ) 

IF^IDMF^Q 1 ! - AND. ISTART. GE.ITON2) MODE=MODE2 
NSAT = 1 

IF (ID .EQ. ' D ' ) NSAT =2 
IF (ID .EQ. 'E') NSAT=3 

C Read the Vae Data File, bin error-free channel one data, 

C convert to Rayleighs, and store m arrays. 

DO 290 NRE AD =1, NREADMAX 
CALL VAERD 

IF (IE .EQ. 'E ') GOTO 260 
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IF 

IF 

IF 

IF 

IF 

IF 


(JFW NE. 10 .AND. JFW .NE. IFW) GOTO 290 
(OADATA {9 ) . LT . ZMIN) GOTO 290 

lf E ^ A ^!i GT 0R Z ™T T .i 9 S ?0P, STOP 

IF aSQll Q EQ.i°OR.ISQl2.EQ.l)GOTO 270 
IF ( ISQ1 1L . EQ . 1 . OR . ISQ12L . EQ . 1 ) GOTO 270 
IF (IATN1.EQ.1) GOTO 280 
GOTO 242 

241 CONTINUE 

IF (ISQ2 ,EQ.l)GOTO 270 
IF (IATN2 .EQ. l)GOTO 280 

242 CONTINUE 
CALL RAYL 
ISKIP=1 

IF (ICE.NE.l) GOTO 290 

CALL SUNVEC (IDATE, ITIME, SUN, GST) 
sth=oadata (3) 

i f^oadata ( 4 ) . 1 1 . 0 . 0 ) sph=3 60 . 0+oadata ( 4 ) 
gstr=gst/57. 29578 
Tmat (1 , 1) =COS (gstr) 

Tmat (l, 2 )=-l- 0 *SIN(gstr) 

Tmat ( 1 , 3 ) =0 . 0 
Tmat (2,1) =SIN{gstr) 

Tmat (2,2) =COS (gstr) 

Tmat (2 , 3 ) =0 . 0 
Tmat (3 , 1 ) =0 . 0 
Tmat ( 3 , 2 ) =0 . 0 
Tmat (3 , 3 ) =1 . 0 
CALL CART(sth,sph,sgeo) 

CALL MAXMULl(Tmat,sgeo,sgei) 
sgei (l)=sgei (1) * (oadata 1) +6370.0) 

svMRayui) UwU2)' (3) +Rayl ( 4 ) ) /4 .0 

IiJSqRt'iI . 0* (ICH11+ICH12+ICH11L+ICH12L) ) *STV(1) /4 .0 
za=ABS ( OADATA ( 9 ) ) 

GOTO 246 

245 s ig =SQRT* ABS ( 1 . 0 * ICH2 ) ) * STV ( 2 ) 

246 ■ 

> 1E13.5,’ ‘ , 1E13 . 5 , 1 , 1E13 .5) 

GOTO 290 

260 NERR = NERR + 1 
GOTO 290 

270 NSQ = NSQ + 1 
GOTO 290 

280 NATN = NATN + 1 
GOTO 290 
290 CONTINUE 
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STOP 

END 


SUBROUTINE VAERD 

C This version reads VAX formatted VAE data. 

C CALLING PROGRAM MUST PROVIDE COMMON /CVAERD/ AND PRESET IEND4 
C TO ZERO . IEND4 IS RESET TO 1 UPON EOF ON VDF . IF CALLING ROUTINE 
C DOES NOT CHECK THIS IT MAY GO INTO AN INFINITE LOOP UPON EOF. 

DIMENSION REC (295) ,IREC(295) ,IBUFF(295) ,BUFF{295) ,OADATA(22) , 

> Al(3) ,A2 (3) ,B1(3) ,B2 (3) ,E(3) ,F(3) ,P(3) ,DF1 (3) ,DF2 (3) , 

> R(3) ,V1 (3) ,V2 (3) ,V1L(3) ,V2L(3) , RMOON(3) , SUN(3) 

LOGICAL* 1 SATID ( 4 ) , ID , SATIDl , SATID2 , SD1 , SD2 , NI 1 , NI 2 
LOGICAL*4 MODE1 ,MODE2 

CHARACTER* 40 INDNAME 
CHARACTER *4 IE 
INTEGER* 4 SHFTR 

INTEGER* 2 LEN _ 

EQUIVALENCE (REC, IREC) , (IBUFF,BUFF) , (SATID, IREC (3 ) ) 

COMMON / CVAERD/ ITIME , THET, ICH11, ICH12 , ISQ11, ISQ12, ITIMEL,THETL, 

> IFW, IE, ICH11L, ICH12L, ISQ11L, ISQ12L, IATN1 , ICH2 , ISQ2 , IATN2 , 

> TBAF1 , TBAF2 , TBAF11 , TBAF22 , TAEL , TBEL , TPMT1 , TPMT2 , TFW , IEND4 , 

> IX, ISVIF, OADATA,OADATL 


COMMON /CVEC/ Vl,V2,VlL,V2L,R,sun,gst 


DATA P/0 . , - .398, .917/ , LORBIT/O/ 

DATA SAVE1/2 . 0/ , PI/3.14159/, FAC/57.29578/ 
DATA ISW, IIREAD/2*0/, INUNIT/4/, IANGSW/0/, 


PI2/6. 28318/ 


C READ HEADER AND FIRST RECORD UNLESS FILE IS PRE POSITIONED USING 
C ENTRY INITV 


IF ( ISW .NE. 0) GO TO 100 

READ ( INUNIT, END=2 50, ERR=2 50) (IREC (I) ,1=1,12) 


IDATE=IREC (1) 
IORBIT=IREC (2) 
ISF=1 


READ ( INUNIT , END=250 , ERR=250) (IBUFF(K) ,K=1,295) 


ISW=1 

100 IIREAD=IIREAD+1 

IF (IIREAD.NE.l) GO TO 400 
DO 120 K=1 ,295 

120 IF^REC (275y F GT K 23 .0 . AND . BUFF (275) . LT . 1 . 0 ) REC (275 ) =REC (275 ) -24 . 0 
IF (REC (277) . GT .179. 0 . AND . BUFF (277 ) .LT . -179 . 0) REC (277) =REC (277 ) - 
>360.0 

IF (REC (277 ) .LT.- „„ „ 

>179 . 0 .AND. BUFF (277) .GT . 179 . 0)REC (277) =REC (277) +360.0 
IF (REC (281) . GT .179.0. AND . BUFF (281) . LT .-179.0) REC (281) =REC (281) - 
>230.0 


READ ( INUNIT , END=250 , ERR=250 ) (IBUFF(K) ,K=1,295) 
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C IF ORBIT IS SPINNING, FIND OA VECTORS NEEDED FOR INTERPOLATIONS 
IEM ISfTnE . 1 . OR . ISVIF . EQ . 0 ) GOTO 2 60 

190 S(Re2(kI 2 l?:Tooo.or.rec(k).gt.iooo,iaflag = o 

IF ( I AFLAG . EQ . 0 ) GOTO 260 

IF (LORBIT . EQ . IORBIT) GOTO 200 

LORBIT=IORBIT 

CALL SUNVEC (IDATE, ITIME, SUN,gst) 

CALL CROSS (P, SUN, E) 

200 CALL CART (REC (287) ,REC(286) ,A1) 

CALL CART (REC ( 293 ), REC (292) ,A2) 

CALL CART (BUFF (287) , BUFF (286) ,B1 
CALL CART (BUFF (293) , BUFF (292 ), B2 ) 

DO 210 NN=1,3 

DF1 (NN) =A1 (NN) -B1 (NN) 

210 DF2 (NN) =A2 (NN) -B2 (NN) 

CALL CROSS (DF2,DF1,R) 

CALL DETER ( A1 , Bl , R , DETl ) 

CALL DETER ( A2 , B2 , R , DET2 ) 

GAMMA=PI2 -ANGLE ( A1 , Bl ) 

CMCHl=COS (REC (285) /FAC) 

CMCH2=COS (REC (291) /FAC) 

CALL LUNVEC ( A1 , A2 , CMCH1 , CMCH2 , RMOON) 

GO TO 260 

C IF EOF, SET FLAG AND RETURN 


250 


IEND4=1 
GO TO 4 


C UNPACK THE DATA RECORD 


260 


270 

280 

290 


300 


. LE . 1.0 ) GO 


IF ( REC (15) 

IF ( SAVE1 .NE. 2 ) 
REC (16) = 1.570681 
REC (17) =0.0 
REC (15) =1.0 
GO TO 290 
REC (15)= SAVE1 
GO TO 300 

REC (15) 


GO TO 


TO 280 
270 


SAVE1 

DELTH 

THET1 

THET2 

ITIME 

ITIMEL 

TBAF1 

TBAF2 

TBAF11 

TBAF22 

TAEL = 

TBEL 

TPMT1 


.125 * REC (17) 
REC (16) - DELTH 

REC (16) + .0626 
IREC(l) - 125 
= I REC { 1 ) - 63 
= REC (3) 

= REC (4) 

= REC (5) 

= REC (6) 

REC (7) 

REC (8) 

= REC (9) 


* REC (17) -DELTH 


TPMT2 = REC (10) 

TFW = REC (11) 

C JJ = INDEX TO STATUS BITS 


400 


JJ = IIREAD *4+14 

ITIME = ITIME + 125 

ITIMEL = ITIMEL + 125 

IF ( IREC ( JJ+1 ) -NE. -1 ) GO TO 410 
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ICH11 = -1 
ICH12 = -1 
GO TO 420 

410 ICH11 = JIBITS (IREC (JJ+1) ,16,16) 

ICH12 = JIBITS(IREC(JJ+1) ,0,16) 

420 IF ( IREC ( JJ+2 ) .NE. -1 ) GO TO 430 
ICH11 = -1 
ICH12 = -1 
GO TO 440 

430 ICH11L = JIBITS (IREC (JJ+2) , 16, 16) 

ICH12L = JIBITS ( IREC (JJ+2) ,0,16) 

440 IF ( IIREAD . NE . 1 ) ICH2 = IREC(JJ+3) 

ITEMP = IREC (JJ) 

C ZENITH ANGLE CALCULATIONS (OLD METHOD) 

THET =0.0 
THETL =0.0 

IF (THET1 . GT . PI2 .OR. THET2.GT.PI2 .OR. REC (15) .GT. 1 . 0)GO TO 450 
IF ( IANGSW .EQ. 1) GO TO 450 
THET1 = THET1 + DELTH 

THET = ACOS (COS (THETl ) * REC (15)) * 57.29583 
IF (THET1 .LT. 0.0 ) THET = 360.0 - THET 
IF (THETl .GE. 6.28318) THETl = THETl - 6.28318 
IF (THETl .GT. 3. 14159) THET = 360.0 - THET 
THET 2 = THET 2 + DELTH 

THETL = ACOS (COS (THET2) * REC(15)) * 57.29583 
IF (THET2 .LT. 0.0) THETL = 360.0 - THETL 
IF(THET2 .GE. 6.28318) THET2 = THET2 - 6.28318 
IF (THET2 .GT. 3.14159) THETL = 360.0 - THETL 
450 CONTINUE 

C UNPACK STATUS WORD 

IWORD= JIBITS ( ITEMP , 3 1 , 1 ) 

IF ( IWORD .EQ. 0 ) GO TO 500 
IE= 1 E ' 

GO TO 501 

500 IE= ' 

501 ISQ2=JIBITS ( ITEMP, 24,1) 

ISQ11= JIBITS { ITEMP, 19,1) 

ISQ12=JIBITS (ITEMP, 18,1) 

ISQ11L=JIBITS (ITEMP, 17,1) 

ISQ12L=JIBITS (ITEMP, 16,1) 

IATN1=JIBITS (ITEMP, 9 , 1) 

IATN2=JIBITS (ITEMP, 8, 1) 

IFW=JIBITS ( ITEMP ,0,3) 

C PERFORM OA INTERPOLATIONS 

RAT= (IIREAD- 1) /64 . 

RATL= (IIREAD- . 5) /64 . 

DO 700 N=1 , 8 

O AD AT A (N) =REC (273+N) + (BUFF (2 73 +N) -REC (273 +N) ) *RAT 
700 OADATL (N) =REC (273+N) + (BUFF (273+N) -REC (273+N) ) *RATL 
IF (O AD AT A ( 2 ) . LT . 0 . 0 ) OADATA ( 2 ) =OADATA ( 2 ) +2 4 . 0 
IF (OADATA ( 2 ) . EQ . 0 . 0 ) OADATA ( 2 ) =2 4 . 0 
IF (OADATL ( 2 ) . LT . 0 . 0 ) OADATL ( 2 ) =OADATL ( 2 ) +24 . 0 
IF (OADATA (8) .LT. -180 . 0) OADATA (8) =OADATA ( 8 ) +360 . 0 
IF (OADATL (8) . LT . -180 . 0) OADATL (8) =OADATL (8) +360 . 0 

IF ( ISF . EQ . 1 ) GOTO 720 
DO 710 N=9 , 22 
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710 

720 


O AD AT A (N) =REC (273+N) + (BUFF (273+N) -REC (273+N) ) *RAT 
OADATL (N) =REC (273+N) + (BUFF (273+N) -REC (273+N) ) *RATL 
GOTO 800 

OADATA ( 9 ) =REC (282 ) +RAT*8 . *REC (283 ) 

IF (OADATA (9 ) .GE . 180 . ) OADATA (9 ) =OADATA (9 ) -3 60 
IF (OADATA (9) . LE . -180 . ) OADATA (9 ) =OADATA (9 ) +360 . 
OADATA (10) =REC (283) 

OADATL (9) =REC (282 ) +RATL*8 . * REC (283) 

IF (OADATL ( 9 ) . GE . 1 8 0 . ) OADATL ( 9 ) =OADATL ( 9 ) -3 6 0 . 

IF (OADATL (9) . LE . -180 . ) OADATL (9 ) =OADATL (9 ) +360 . 
OADATL (10) =REC (283) 

IF ( ISVIF . EQ . 0 ) GOTO 800 
IF ( I AFLAG . EQ . 0 ) GOTO 750 
DELTA=RAT *GAMMA 
CD=COS (DELTA) 

CE=COS (GAMMA-DELTA) 

CALL SOLVE (A1 , B1 , R, CD, CE, 0 , DET1 , VI ) 

CALL SPHERE (VI, OADATA (14) , OADATA (13) ) 

OADATA (11) = ANGLE (VI , SUN) *FAC 
OADATA (12) = ANGLE (VI , RMOON) *FAC 
CALL MULT (SUN, E, P, VI , F) 

CALL S PHERE ( F , OADATA (15), OADATA (16)) 

CALL SOLVE (A2 , B2 , R , CD, CE , 0 , DET2 , V2 ) 

CALL S PHERE (V2 , OADATA (20), OADATA (19)) 

OADATA ( 17 ) = ANGLE (V2 , SUN) *FAC 
OADATA ( 18 ) = ANGLE (V2 , RMOON) *FAC 
CALL MULT (SUN, E,P,V2,F) 

CALL SPHERE (F, OADATA (21) ,OADATA{22)) 


DELTAL=RATL*GAMMA 
CDL=COS (DELTAL) 

CEL=COS (GAMMA -DELTAL) 

CALL SOLVE (A1 , B1 , R,CDL, CEL, 0 , DET1 , VI L) 
CALL SPHERE ( VI L, OADATL (14) , OADATL ( 13 ) ) 
OADATL (11) = ANGLE ( V1L , SUN) *FAC 
OADATL (12) = ANGLE (V1L, RMOON) *FAC 
CALL MULT (SUN, E,P,V1L,F) 

CALL SPHERE (F, OADATL ( 15 ), OADATL (16) ) 
CALL SOLVE (A2 , B2 , R, CDL, CEL, 0 , DET2 , V2L) 
CALL SPHERE (V2L, OADATL (20) , OADATL (19) ) 
OADATL ( 17 ) = ANGLE (V2L, SUN) *FAC 
OADATL (18) = ANGLE (V2L, RMOON) *FAC 
CALL MULT (SUN, E,P,V2L,F) 

CALL SPHERE ( F , OADATL (21), OADATL (22)) 
GOTO 800 


750 DO 755 N=ll,22 

OADATA (N) =-99999 . 

755 OADATL (N) =-99999 . 

C APPLY DEAD TIME CORRECTION 

800 IF ( ICH11 .LT. 0) GO TO 810 

ICH11=ICH11/ (1-ICHll/l . 11E5) 
ICH12=ICH12/ (1-ICH12/1 . 11E5) 

810 IF (ICH11L .LT. 0) GO TO 820 

ICH11L=ICH11L/ (1-ICHllL/l . 11E5) 
ICH12L=ICH12L/ (1-ICH12L/1 . 11E5) 
820 IF ( ICH2 .LT. 0) GO TO 830 
ICH2=ICH2 / (1-ICH2/4.44E5) 

830 IF (IIREAD.EQ. 64) IIREAD=0 
IX=IIREAD 
4 RETURN 
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179 

900 

910 


920 


940 

950 


ENTRY INITV (IUNTP, JTIMEP, IDDATE, IORB, ID, ISPINF, INF, 
>MODEl , ITON, ITOFF, IDMF,MODE2 , ITON2 , ITOFF2) 


INUNIT= IUNTP 

I I READ = 0 

ISW = 0 

SAVE1 =2.0 

JTIME= JTIMEP - 8000 

READ ( INUNIT, END=4, ERR=4) (IREC(I) ,1=1,12) 
ID = SAT ID ( 1 ) 

IDDATE = IREC(l) 


IDATE= IDDATE 
IORB = IREC (2 ) 

IORBIT=IORB 

IF ( ID . EQ . ' C 1 . OR . ID . EQ . ' c * ) INDNAME= 

+ 1 SPRLC$DISK1 : [VAECOMMON . IND] VDFC . DAT ' 

IF ( ID . EQ . ' D ' . OR . ID . EQ . ' d ' ) INDNAME= 

+ ' SPRLC$DISK1 : [VAECOMMON .IND] VDFD . DAT 1 
IF ( ID . EQ . ' E ' . OR . ID . EQ . ' e 1 ) INDNAME= 

+ ' SPRLC$DISK1 : [VAECOMMON . IND] VDFE . DAT ' 

OPEN (UNIT=9 , NAME=INDNAME, STATUS= 1 OLD' , READONLY , ERR=179 ) 
GOTO 900 

TYPE * , ' ERROR OPEN 1 

READ (9,910, END=9 40 ) SATID1 , IORB1 , MODE1 , SD1 , Nil , IDATE1 , 
+ITON1 , ITOFF1 , MT1 

FORMAT (Al , IX, 15 , IX, A4 , IX, 2A1 , IX, 3 ( 15 , IX) , 13 ) 

IF ( IORB1 . NE . IORB) GOTO 900 
IORB2=0 

READ (9,910, END=92 0 ) SATID2 , IORB2 , MODE2 , SD2 , NI2 , IDATE2 , 
+ITON2 , ITOFF2 , MT2 


ISF=0 

IF (SD1.EQ. 'S' ) ISF=1 

ISPINF=ISF 

INF=0 

IF (Nil . EQ . 'N' ) INF=1 

ITON=ITONl 

ITOFF=ITOFFl 

IDMF=0 

IF { IORB1 . EQ . IORB2 ) THEN 

IDMF=1 

ELSE 

MODE 2= 1 

ITON2=0 

ITOFF2=0 

ENDIF 

CLOSE (UNIT=9) 

GOTO 950 

TYPE *, 'ORBIT NOT FOUND IN INDEX' 

CLOSE (UNIT=9) 

STOP 

CONTINUE 

READ ( INUNIT, END=4,ERR=4) (IBUFF(K) ,K=1,295) 
ISW= ISW+1 

IF (IBUFF ( 1 ) .LT. JTIME) GO TO 950 
RETURN 

ENTRY SETUP ( IANGZ) 

IANGSW= IANGZ 
GO TO 4 
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END 


FUNCTION ANGLE (VI, V2) 

ftNrfF-ACOs'nvlll^*V2 (1) +V1 {2 ) *V2 (2) +V1 (3 1 *V2 (3 ) ) / 

>SQRT ( (VI ( 1 ) * *2+Vl (2 ) **2+Vl (3) **2) * (V2 (1) **2+V2 (2 ) **2+V2 (3 ) * *2 ) ) ) 

RETURN 

END 


SUBROUTINE CROSS(A,B,C) 
DIMENSION A(1),B(1),C(1) 

C ( 1 ) =A (2 ) *B (3 ) -A (3 ) *B (2 ) 
C(2)=A(3) *B(1) -A(l) *B(3) 

C (3 ) =A (1) *B(2) -A (2) *B(1) 

CL=SQRT (C(l) **2+C(2) **2+C(3) 

C(1)=C(1) /CL 

C(2)=C(2) /CL 

C(3)=C(3) /CL 

RETURN 

END 


** 2 ) 


SUBROUTINE MULT (SUN, E,P,V,F) 

DIMENSION SUN ( 1 ) , E ( 1 ) , p ( 1) , v ( 1 ) < F(1 ) 
p ( i) =SUN { 1 ) *V { 1 ) +SUN ( 2 ) *V (2 ) +SUN (3 ) *V (3 ) 
F (2 ) =E ( 1) *V ( 1 } +E (2 ) *V (2 ) +E (3 } *V (3 ) 

F (3 ) = P ( 1 ) *V ( 1 ) + P (2 ) *V (2) +P(3) *V (3 ) 

RETURN 

END 


SUBROUT INE CART ( THETA , PH I , W ) 
DIMENSION W(l) 

T=THETA* .01745329 
P=PHI* . 01745329 
W(l) =COS (T) *COS (P) 

W (2 } =COS (T) *SIN (P) 

W (3 ) =SIN (T) 

RETURN 

END 


SUBROUTINE SPHERE (V, TH, PH) 
DIMENSION V{1) 

AV=SQRT (V(l)**2+V(2)**2+V(3)**2) 
V ( 1 ) =V ( 1 ) / AV 

V ( 2 ) = V ( 2 ) / AV 

V ( 3 ) = V ( 3 ) / AV 

PH=ATAN2 (V{2) ,V(1) ) *57 .29578 
TH=ASIN (V { 3 ) ) *57 .29578 
RETURN 
END 


100 


SUBROUT INE SOLVE ( X , Y , Z , CAX , C AY , CAZ , DET , V) 

DIMENSION X(1),Y(1),Z(1),V(1) 

•(C&*-(? , (2?*?(3?-Y(3)*Z(2)) ♦ X ( 2 ) * (Y (3 ) *CAZ-CAY*Z (3 ) 

V{2)- ] UalMCAY*Z(3)-?^3)*CAzH CAX* (Y (3 ) *Z ( 1) -Y ( 1) *Z (3 ) 

X( 2 )*(C A V*Z(1>- Y (1)*C AZ , 

+ CAX* (Y(1)*Z(2)-Y(2)*Z{1))) / DET 


) 

) 


RETURN 
V(1)=0. 
V ( 2 ) =0 . 
V (3 ) =0 . 
RETURN 
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END 


SUBROUTINE DETER (X, Y, Z, DET) 
DIMENSION X(l) ,Y(1) ,Z(1) 

DET= X(1)*{Y(2)*Z(3)-Y(3)*Z(2)) 
> + X(3)*(Y(1)*Z(2)-Y(2)*Z(1)) 
RETURN 
END 


+ X(2)*(Y(3)*Z(1)-Y(1)*Z(3)) 


SUBROUTINE LUNVEC (A1 , A2 , G, H, RMOON) 

DIMENSION A1 ( 1 ) , A2 ( 1 ) , RMOON ( 1 > 

A=A1 (1) 

B=A1 (2) 

C=Al (3) 

D=A2 (1) 

E=A2 (2) 

F=A2 (3) 

S = (E*C-B*F) **2 + (A*F-D*C) **2 + (E*A-D*B)**2 
q = (A*H-D*G) * (F*A-D*C) + (E*G-B*H) * (E*C-B*F) 

T = (E*G-B*H) **2 + (A*H-D*G) **2 - (E*A-D*B) **2 

R = (2*Q) **2 - 4*S*T 

IF ( R . LT . 0 . ) GOTO 100 

RMOON ( 3 ) = (2*Q+SQRT (R) ) / <2*S) 

V = H-F*RMOON (3 ) 

W = G-C*RMOON (3 } 

RMOON (1) = (W*E-V*B) / (E*A-D*B) 

RMOON (2) = (V*A-D*W) / (E*A-D*B) 

RETURN 

100 RMOON ( 1 ) = 0 . 

RMOON (2) =0. 

RMOON (3) =0. 

RETURN 

END 


SUBROUTINE SUNVEC (IDATE, ITIME , SUN, gst ) 

DIMENSION SUN { 1 ) 

DATA RAD/57.29578/ 

REAL *8 DJ , FDAY 
FDAY=FLOAT (ITIME) /86400000 . 

IYR=IDATE/1000 

IDAY=IDATE-IYR*1000 

DJ=365*IYR+ (IYR-1) / 4+IDAY+FDAY-O . 5D0 
T=DJ/36525 . 

VL=DMOD(279 . 696678+. 98564733 54*DJ, 3 60. DO) 

gst=dmod (279. 690983 +. 98564733 54*DJ+360.*FDAY+180. ,3 60. DO) 

G=DMOD (3 58 .475845+. 985600267*DJ, 360. DO) /RAD 

SLONG=VL+ (1 .91946- ,004789*T) *SIN (G) + . 020094*SIN (2 . *G) 

OBLIQ= (23 .4 52 29-0. 0130125 *T) /RAD 

SLP= (SLONG- . 005686) /RAD 

SIND=SIN (OBLIQ) *SIN (SLP) 

COSD=SQRT ( 1 . -SIND* *2 ) 

SDEC=ATAN(SIND/COSD) 

SRASN=3 . 14159-ATAN2 ( 1 /TAN (OBLIQ) *SIND/COSD, -COS (SLP) /COSD) 
SUN ( 1 ) =COS ( SRASN ) *COS (SDEC) 

SUN (2 ) =SIN ( SRASN) *COS (SDEC) 

SUN ( 3 ) =S IN ( SDEC ) 

RETURN 

END 


SUBROUTINE RAYL 

C SUBROUTINE RAYL CONVERTS RAW VAE COUNTS TO RAYLEIGHS, AND STORES 
C THE VALUES IN RAY1(1:4) AND RAY2 . THE CHANNEL NUMBER REQUESTED 
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C MUST BE IN ICH AND THE SATELLITE NUMBER {C=l , D=2 , E=3 ) MUST BE 
£ first DARK COUNT IS ESTIMATED FROM THE 

C sSuS'aS) SSct/SJt' 'S fSctlC AND ZODIACAL BACKGROUND 
c IF IGZC=1 . 

C CALLING PROGRAM MUST SUPPLY /CVAERD/ AND /CRAYL/ IN COMMON 

C Modification 6/15/88 to work with galactic and zodiacal background 

c qubtract ion on SPRLC. . , * 

c * errors detected in VAEREAD on 1/7/91 fixed here * 

C Caning program must initialize ISKIP to 0 and then set it to 1 
C immediately after RAYL is called for the first time. 

INTEGER* 2 LEN 

INTEGER* 4 RA ( 4 ) , DEC ( 4 ) , ELAT ( 4 ) , ELON { 4 ) 

CHARACTER* 4 IE A „, c „ 

DIMENSION GAL (5, 120, 60) ,ZOD (5,60,30), 

> OA(22) ,OAL(22) ,RAY1(4) ,IFM(7,2,3) ,FBW(7,2,3) 

nTMFKTSTDN Cl (2 3) C2 (2,3) ,TPM(2) ,OLDTPM(2) ,DC(2) ,s(7, 7,2,3) , 

^DIMENSIO^Cl (2 . 3 > ; « < l ; ^ , 2 ; , OLDTFW ( 2 ) . IOLDFW ( 2 , 

EQUIVALENCE (C, S (1) ) , (D, S (99) ) , (E,S(197)) 

COMMON^ /CVAERD/ ITIME^THET^ICHll^ICHl 2^ ISQll^ISQ^^ITIMELiTHETL; 

TBAF1 , TBAF2 , TBAF11 , TBAF22 , TAEL, TBEL, TPM, TFW 
, IEND4 , IX , ISVIF , OA, OAL 


> 

> 

> 


COMMON /CRAYL/ ICH , NSAT , IGZC , STV , DC , RAY1 , RAY 2 , ICE , ISKIP 

DATA OLDTPM/0 . , 0 . / , OLDTFW/0 . , 0 • / , IOLDEW/0 0/ 

DATA Cl/. 1544, .175, .1658, .1677 ,*1757 ‘ / 

DATA C2/4. 5,4. 6,19. 02 ,11.85,11.57, -1.827/ 


2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 
7 

D. 

2 

3 

4 

5 

6 
7 
1 
2 

3 

4 

5 

6 


C / 

20. C 

1 . 

19.6 

> r 

19.4 


19.4 

^ / 

19. S 

>, 

21.8, 24.4, 

101 . 5 . 

79.3, 

66.8, 

57.3, 

51.1, 

, 45.9, 

, 44 . 2 . , 

12 

. 2 , 

11 

.5, 

11 

• 1, 

11 

.0, 

11 

.1, 

12 

.0, 

13.5, 

23 

. 6, 

20 

.7, 

19 

.2, 

18 

.6, 

18 

.4, 

18 

. 6 , 

19.6, 

35 

. 2 , 

35 

.3 , 

35 

.8, 

36 

.7, 

37 

.9, 

39 

.4, 

41.4, 

0 . 

0 , 

0. 

0, 

0. 

0, 

0. 

0, 

0. 

0, 

0. 

0, 

0.0, 

10 

.7 , 

10 

.4, 

10 

.6, 

10 

.9, 

11 

• 4, 

12 

• 1, 

12.8, 

0 . 

30, 

0. 

27, 

0. 

25, 

0. 

24, 

0. 

24, 

0. 

24, 

0.25, 

0 . 

00, 

0 . 

00, 

0. 

00, 

0. 

00, 

0. 

00, 

0. 

00, 

0.00, 

1 . 

45, 

1 . 

13, 

0. 

95, 

0. 

82, 

0. 

73, 

0. 

65, 

0.63, 

0 . 

15, 

0. 

15, 

0. 

15, 

0. 

16, 

0. 

17, 

0. 

17, 

0.19, 

0 . 

24 , 

0. 

24, 

0. 

23, 

0. 

23, 

0. 

24, 

0. 

26, 

0.29, 

0 . 

15, 

0 . 

14, 

0. 

14, 

0. 

14, 

0, 

14, 

0. 

15, 

0.17, 

0 . 

61 , 

0 . 

61, 

0 . 

62, 

0. 

64, 

0. 

66, 

0. 

68, 

0.72 / 

D/ 

80 . i 

3/ 

72 . 

1, 

68. 1 

5, 

72. 

1, 

77 . 

8, 

87 . 

3,106.6 

7 . 

6, 

7 . 

6, 

7 . 

9, 

8. 

4, 

9. 

6, 

12. 

2, 

16.3, 

26 

5.0, 

24 

1.2, 

23 

5.3, 

23 

1.5, 

24 

1.4, 

26.2, 

29.2, 

9 . 

7 , 

8 . 

5, 

8 . 

2, 

8. 

1, 

8. 

1, 

8. 

1, 

8.3, 

12 

> .2, 

11 

..3, 

11 

L . 1, 

11.1, 

n.i. 

11.2, 

11.6, 

n 

O . 

0 . 

.0. 

0, 

, 0, 

0. 

0, 

0. 

0, 

0. 

.0, 

0.0, 

2 2.2. 

17.6, 

16.7, 

16.3, 

16.4, 

19.7, 

20.4, 

0 , 

. 13 , 

0 , 

.12, 

0 . 

.11, 

0. 

.11, 

0. 

.11, 

0, 

.11, 

0.11, 

0 . 

.00, 

0, 

.00, 

0. 

.00, 

0. 

.00, 

0, 

.00, 

0, 

.00, 

0.00, 

0 

. 10, 

0 

.10, 

0 

.10, 

o. 

.11, 

0, 

.12, 

0. 

.16, 

0.21, 

0 

.24, 

0 

.19, 

0 

.18, 

0 

.18, 

0 

.18, 

0 

.21, 

0.22, 

0 

.76, 

0 

.69, 

0 

.65, 

0 

.69, 

0 

.74, 

0 

.83, 

1.02, 

0 

.37, 

0 

.34, 

0 

.33, 

0 

.33, 

0 

.34, 

0 

.37, 

0.41, 
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7 

0 . 

14 , 

0 . 

13, 

0 . 

13, 

DATA 

E/124.4 

, 

88.8 

t 

73. 

2 

11 

. 0 , 

10 

.3, 

10 

.1, 

3 

42 

.9, 

42 

• 9, 

42 

.8, 

4 

22 

.8, 

21 

.2, 

20 

.4, 

5 

12 

.7, 

12 

.7, 

13 

• 2, 

6 

0 . 

0 , 

0 . 

0 , 

0 . 

0 , 

7 

19 

.8, 

19 

• 6, 

20 

.4, 

1 

0 . 

26, 

0 . 

24, 

0 . 

23, 

2 

0 . 

00, 

0 . 

00, 

0 . 

00, 

3 

0 . 

15, 

0 . 

14, 

0 . 

14, 

4 

0 . 

23, 

0 . 

23, 

0 . 

24, 

5 

1 . 

27, 

0 . 

91, 

0 . 

76, 

6 

0 . 

60, 

0 . 

60, 

0 . 

60, 

7 

0 . 

16, 

0 . 

16, 

0 . 

16, 


0 . 

13, 

0 . 

13, 

0 . 

13, 

0 . 

13 / 


61.7 

I 

53.9 

f 

47 . 

2 , 

44.2, 

io 

• 1, 

10 

.4, 

11 

.7, 

14 

.1, 

42 

• 9, 

42 

.9, 

42 

.9, 

42 

.9, 

20 

.7, 

21 

.3, 

22 

.8, 

29 

.3, 

15 

.8, 

21 

.8, 

29 

• 3, 

42 

• 4, 

0 . 

0 , 

0 . 

0 , 

0 . 

0 , 

0 . 

0 , 

22 

.5, 

24 

• 2, 

46 

.2, 

72 

• 8, 

0 . 

23, 

0 . 

24, 

0 . 

26, 

0 . 

33, 

0 . 

00, 

0 . 

00, 

0 . 

00, 

0 . 

00, 

0 . 

14, 

0 . 

14, 

0 . 

16, 

0 . 

19, 

0 . 

26, 

0 . 

29, 

0 . 

54, 

0 . 

86, 

0 . 

63, 

0 . 

55, 

0 . 

48, 

0 . 

28, 

0 . 

60, 

0 . 

60, 

0 . 

60, 

0 . 

60, 

0 . 

19, 

0 . 

27, 

0 . 

36, 

0 . 

52 / 


DATA 

IFM/ 

55, 73, 52, 63, 33, 0, 42, 


> 

63, 

0, 73, 42, 55, 52, 33, 



> 

73, 

52, 42, 48, 55, 0, 63 

/ 


> 

48, 

0, 52, 63, 73, 42, 55, 



> 

73, 

52, 28, 65, 55, 0, 63 

, 


> 

65, 

0, 52, 63, 73, 28, 55 / 



DATA 

FBW/ 

29.8, 20.0, 22.5, 22.7, 

0,0, 

o.o, : 

> 

22 . 

7, 0.0, 20.0, 23.0, 29.8, 

22.5, 

0.0, 

> 

29 . 

0, 21.0, 19.0, 0.0, 25.0, 

0.0, 

20.5, 

> 

0.0 

, 0.0, 21.0, 20.5, 29.0, 

19.0, 

25.0, 

> 

15 . 

8, 20.1, 0.0, 0.0, 19.5, 

0.0, 

21.0, 

> 

0.0 

, 0.0, 20.1, 21.0, 15.8, 

0.0, 

19.5 / 


100 


ICE=0 

IF (IFW. LT . 1 .OR . IFW.GT . 7 )GOTO 300 
IF(TPM(ICH) . EQ . OLDTPM ( ICH) ) GOTO 100 

IF (TPM ( ICH) . GT . 40 . 0 . OR . TPM ( ICH) . LT . -2 0 . 0 ) TPM ( ICH) = OLDTPM ( ICH) 

OLDTPM ( ICH) =TPM( ICH) . 

DC (ICH) =EXP (Cl (ICH,NSAT) * (TPM (ICH) -C2 (ICH,NSAT) ) ) 


IF (TFW. EQ .OLDTFW (ICH) . AND . IFW . EQ . IOLDFW ( ICH) )GOTO 200 
IF (TFW .GT .40.0. OR . TFW . LT . -20 . 0) TFW=OLDTFW ( ICH) 

OLDTFW (ICH) =TFW 
IOLDFW (ICH) = IFW 

S1=S ( IFIX (TFW+3 0 . ) /10 , IFW, ICH,NSAT) 

S2=S ( IFIX (TFW+40 . ) /10, IFW, ICH,NSAT) 
STV(ICH)=S1+(TFW/10. -IFIX (TFW) /10) * (S2-S1) 


200 IF ( ICH . EQ . 2 ) GOTO 250 

RAY1 (1) = (ICH11-DC (1) ) *STV ( 1 ) 
RAY1 (2 ) = ( ICH12 -DC ( 1 ) )*STV(1) 
RAY1 (3) = (ICH11L-DC (1) ) *STV(1) 
RAY1 (4) = (ICH12L-DC (1) ) *STV(1) 
GOTO 500 


250 RAY2= (ICH2-DC (2) ) *STV (2) 
GOTO 500 


300 RAY2=0 

DO 400 1=1,4 
400 RAY1 (I) =0 . 

RETURN 

500 IF (IGZC.EQ.0) RETURN 

C Only open and read the data file the first time RAYL is called. 
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IF (ISKIP . EQ . 0) CALL GZREAD (GAL , ZOD) 
IF { ICH . EQ . 2 ) GOTO 600 
LAM=IFM ( IFW, ICH, NSAT) 

IF (LAM. NE . 0) GOTO 530 

ICE=4 

RETURN 

530 CONTINUE 


ICE=1 

IF (OA ( 13 ) . LT .-180. .OR.OA(13) .GT . 180 . ) ICE=2 
IF (OA ( 14 ) . LT . -90 . .OR .OA ( 14 ) .GT . 90 . ) ICE=2 
IF (OA ( 15 ) . LT . - 9 0 . . OR . OA ( 1 5 ) .GT.90.) ICE=2 
IF (OA ( 16 ) . LT . -180 . . OR . OA (16) .GT . 180 . ) ICE=2 
IF ( ICE . EQ . 2 ) RETURN 
IF (OA ( 13 ) . LT . 0 . ) OA ( 13 ) =OA ( 13 ) +3 60 . 

IF (OAL (13 ) . LT . 0 . ) OAL (13 ) =OAL (13 ) +3 60 . 

IF (OAL ( 13 ) . GT .350. ,AND.OA(13) . LT . 10 . ) OA ( 13 ) =OA (13 ) +360 . 
IF (OAL (13) . LT .10. .AND.OA(13) .GT.350. ) OAL ( 13 ) =OAL (13 ) +360 
DRA=OAL ( 13 ) -OA ( 13 ) 

DDEC=OAL ( 14 ) -OA(14) 

DELAT=OAL (15) -OA(15) 

OAL (16) =ABS (OAL (16) ) 

OA{16) =ABS (OA (16) ) 

DELON=OAL ( 1 6 ) -OA ( 1 6 ) 


C * the following is changed in this program 1/7/91 
RA (1 ) =IFIX (OA ( 13 ) +0 . 25*DRA) /3 +1 
RA (2 ) =IFIX (OA ( 13 ) +0 . 25*DRA*3 . ) /3 +1 
RA (3 ) = IFIX (OAL ( 13 ) +0 . 25*DRA) /3 +1 
RA ( 4 ) =IFIX (OAL ( 13 ) +0 . 25*DRA*3 . ) /3 +1 
DEC (1) =IFIX (OA (14) +0 . 25*Ddec+90 . ) /3 +1 
DEC (2 ) =IFIX (OA ( 14 ) +0 . 25*Ddec*3 +90.)/3 +1 
DEC ( 3 ) =IFIX (OAL ( 14 ) +0 ,25*Ddec+90 . ) /3 +1 
DEC (4) =IFIX (OAL ( 14 ) +0 . 25*Ddec*3 . + 90.)/3 +1 
ELAT ( 1 ) =IFIX { ABS (OA ( 15) +0 . 25*DELAT) ) /3 + 1 
ELAT (2 ) =IFIX ( ABS (OA ( 15) +0 . 25*DELAT*3 ) ) /3 +1 
ELAT (3 ) =IFIX (ABS (OAL ( 15 ) +0 . 25*DELAT) ) /3 +1 
ELAT (4) =IFIX (ABS (OAL (15) +0 . 25*DELAT*3 ) )/3 +1 
ELON ( 1 ) =IFIX (OA (16 ) +0 . 25*DELON) /3 +1 
ELON ( 2 ) =IFIX (OA ( 16 ) +0 . 25*DELON*3 . ) /3 +1 
ELON (3 ) =IFIX (OAL (16) +0.25*DELON) /3 +1 
ELON (4) =IFIX (OAL (16) +0 . 25*DELON*3 . ) /3 +1 

IF (LAM. EQ . 42 .OR . LAN. EQ . 33 .OR .LAM. EQ . 28) LI = 1 

IF (LAM . EQ . 52 . OR . LAM . EQ .48) LI = 2 
IF (LAM . EQ . 55 ) LI = 3 

IF (LAM. EQ. 63 .OR. LAM. EQ. 65) LI = 4 

IF (LAM.EQ . 73 ) LI = 5 


DO 550 1=1,4 

IF (RA ( I ) . GT .120. OR . RA ( I ) . LT . 1 ) RA ( I ) =1 
IF (ELON ( I ) . GT . 60) ELON ( I ) = 6 0 
IF (ELON ( I ) . LT . 1 ) ELON ( I ) =1 
IF (ELAT (I) .GT. 30) ELAT (I) =30 
IF (DEC ( I ) .GT . 60) DEC (I) =60 

C * following elon(I) & elat(I) replace faulty RA(I) & DEC(I) 
IF (GAL (LI , RA ( I ) ,DEC(I) ) .GT.O.OR. 

>ZOD (LI , elon ( I ) , elat (I) ) .GT.0. ) THEN 
GOTO 545 
END IF 
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c * 
c * 

545 


GAL (LI, RA (I) , DEC ( I ) )=0. 

ZOD (LI , RA ( I ) , DEC ( I ) ) =0 . 

ICE=-1 

original “DEC ( I ) ) *16 . ■ in 545 changed 
oriainal 11 IL" in cont 545 changed to 

RAY 1(1) = RAYl(I) - ( (GAL (LI , RA (I) , DEC^IH ) 

& ZOD { 1 i , ELON ( I ) / ELAT { I ) ) 


i *FBW{IFW, ICH,NSAT) 


550 CONTINUE 
RETURN 


600 


640 


C * 
650 


900 


I +360. 


IF^OAL (19) . LT . - 1 8 0 . . OR . OAL (19) .GT . 180 . ) ICE=2 
IF (OAL (20) . LT . -90 . . OR . OAL (20) .GT . 90 . ) ICE=2 
IF (OAL (21 ) . LT . -90 . . OR .OAL (21) .GT . 90 . ) ICE-2 
IF (OAL (22) . LT .-180. .OR. OAL (22) .GT . 180 . ) ICE=2 
IF ( ICE . EQ . 2 ) RETURN 
IF (OAL ( 19 ) . LT . 0 . ) OAL ( 19 ) =OAL (19 
IRA2=IFIX (OAL (19) ) / 3 +1 
IF (IRA2 .GT . 120) IRA2=120 

IDEC2=IFIX (OAL (20) +90 . ) /3 +1 

IF (IDEC2 .GT . 60) IDEC2=60 
IELAT2=IFIX (ABS (OAL (21) ) ) /3+1 
IF ( I ELAT 2 . GT . 3 0 ) IELAT2 =3 0 
IELON2=IFIX (ABS (OAL (22 ) ) ) /3 +1 
IF (IELON2 . GT . 60) IELON2=60 
LAM=IFM(IFW, ICH,NSAT) 

IF (LAM. EQ. 4 2 . OR . LAM . EQ . 33 . OR . LAM . EQ . 28 ) 


LI = 1 


IF (LAM.EQ . 52 . OR . LAM.EQ .48) 
IF (LAM. EQ. 55) LI = 3 

IF (LAM.EQ. 63 . OR . LAM . EQ . 65 ) 


LI = 2 


LI = 4 


. EQ .73) 

. NE . 0 ) GOTO 


LI 

640 


= 5 


. GT . 0 . )GOTO 650 


IL 


IF (LAM. 

IF (LAM. 

ICE=4 
RETURN 

IF (GAL (LI, IRA2, IDEC2] 

GAL (LI , IRA2 , IDEC2 ) =0 . 

ICE=-1 . . , . .. 

following ” li“ replaces original faulty 

RAY2 = RAY2 - (GAL (LI , IRA2 , IDEC2 ) + ......mx 

& ZOD ( 1 i , IELON2 , I ELAT 2 ) ) *FBW ( IFW , ICH , NS AT) 

RETURN 
ICE=3 
RETURN 
END 

SUBROUTINE GZ READ (GAL , ZOD) 

C Subroutine GZREAD returns a 120x60 array of galaxy data and a 60x30 
C Kra? of zodiacal light data for each of 5 wavelengths. 

S, READONLY) 

READ (12,1000, END=5 ) GZDAT 
1000 FORMAT ( IX, 15F5 . 2 ) 

C There°are 5 pairs of galaxy and zodiacal light data (one for each 

c wavelength) . 

DO 10 I = 1, 5 

C Read galaxy data for one wavelength 

DO 20 K = 1, 60 

DO 30 J = 1, 120 
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IC = IC + 1 

GAL ( I , J , K ) = GZDAT { IC) 

30 CONTINUE 

20 CONTINUE 

C Read zodiacal light data for one 

DO 40 N = 1, 30 
DO 50 M = 1, 60 
IC = IC + 1 
ZOD ( I , M, N) = GZDAT (IC) 

50 CONTINUE 

40 CONTINUE 

10 CONTINUE 

RETURN 
END 

SUBROUTINE MAXMUL1 (A, B,C) 

C IMPLICIT DOUBLE PRECISION (A-H,0 

DIMENSION A (3, 3) ,B(3) ,C(3) 

C ( 1 ) =A ( 1 , 1 ) *B ( 1 ) +A { 1 , 2 ) *B (2 ) +A ( 1 
C (2 ) =A (2 , 1 ) *B ( 1 ) + A (2 , 2 ) *B (2 ) + A (2 
C (3) =A(3 , 1) *B(1) +A(3 , 2) *B (2) +A(3 
RETURN 
END 


wavelength 


-Z) 

, 3 ) *B (3 ) 
,3) *B(3) 
, 3 ) *B (3 ) 
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Appendix 3 

The following is a listing of the Pascal (VAX-11) program used to invert the AE-E 
Visible Airglow Experiment measurements of the 0*C®-2P) twilight airglow emtsston 

at 7320 A. 


program TWI FITTER (input , output); 


{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 

{ 


program for inverting multi-directional twilight 
AE P VAE 7320 A observations developed for the project . 

An Assessment of Twilight Airglow Inversion Procedures 
An Assessmen Atmosphere Explorer observations 

by I.C.McDade 

under NASA Grant NAG 5-1502 

this self-contained VAX-U PASCAL program will return: 


the atomic oxygen scale height, 
the 0*{ 2 P) ionization frequency, 
and the O-atom density @ 250 km, 


H (km) 

Iinf (sec -1 ) 


0250 (cm 3 ) 


that best fit the input 7320-30 A column emission rates 
many of the functions and procedures are taken from: 

Y Cambridge University Press, New York, 1986 

other aspects are based on/he formulation described by 

Me Dade et al . J . Geophys . Res . , Vol 9b, pp. 


} 

} 

} 

} 

} 

} 

} 

> 

} 

} 

} 

} 

} 

} 

I 

} 

} 

} 

} 

} 

} 

} 

} 

} 

} 

} 

} 

> 


function gammln (xx: real) : real; 

const 

stp = 2.50662827465; 
half =0.5; 
one = 1.0; 
fpf = 5.5; 

vax 

x,tmp,ser: double; 
i: integer; 

cof; array [1.. 6] of double; 

begin 

cof [1] := 76.18009173; 
cof [ 2 ] : = -86.50532033; 
cof [3] : = 24 . 01409822 ; 
cof [4] := -1.231739516; 
co f [ 5 ] : = 0 . 120858003e-2 ; 
cof [6] := — 0 . 5 3 63 8 2 e - 5 ; 
x : = XX-one; 
tmp : = x+fpf ; 
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tmp : = (x+half ) * In (trap) -trap; 
ser:= one; 
for j : = 1 to 6 do 
begin 
x : = x+one; 
ser:= ser+cof [j ] /x 

end ; ^ , 

gammln:= sngl (tmp+ln (stp ser) ) 

end; 


procedure gser (a,x: real; var 
label 1; 

const 

itmax = 100; 
eps = 3 . Oe-7 ; 

var 

n: integer; 
sum, del, ap: real; 

begin 

gin : = gammln (a) ; 
if (x <= 0.0) then 

begin 

if (x < 0.0) then 

begin 

wr iteln ( 1 pause in GSER 
readln 


gamser,gln: real); 


x less than O'); 


end; 

gamser : = 0.0 
end 
else 
begin 

ap:= a; 
sum : = 1.0/a; 
del : = sum; 
for n : = 1 to 


itmax do 


begin 

ap : = ap+1 .0; 
del : = del*x/ap; 
sum: = sum+del ; 

if (abs (del ) < abs (sum) *eps) then 
goto 1 

writeln( 'pause in GSER-a too large, itmax too small 1 ) 
readln; 


gamser : = sum*exp {-x+a*ln (x) -gin) 

end 

end; 


procedure gcf <a,X: real; var gammcf,gln: real) 
label 
1 ; 

const 

i tmax = 100; 
eps = 3 . 0e-7 ; 

var 

n: integer; „ _ . 

gold, g , tac , bl , bO , anf , ana, an, a 1 , aO : real; 

begin 

gin: = gammln (a) ; 
gold:= 0.0; 


2 



Final Report NASA Grant NAG 5- 1502 APPENDIX 3 


aO : = 1.0; 
al : = x; 
bO : = 0.0; 
bl:= 1.0; 
fac:= 1.0; 

for n:= 1 to itmax do 

begin 

an:= 1.0*n; 
ana:= an-a; 
a 0 : = (al+a0*ana) *fac; 
b0:= (bl+b0*ana) *fac; 
anf:= an*fac; 
al:= x*a0+anf*al; 
bl:= x*b0+anf*bl; 
if (al <> 0.0) then 
begin 

fac : = 1.0/al; 
g : = bl*fac; 

If (abs ( (g-gold) /g ) < eps) then 
goto 1 ; 
gold:= g 

end 

writeln( ’pause in GCF-a too large, itmax too small’) 
readln; 

gammcf := exp ( -x+a*ln (x) -gin) *g 

end; 


function gammp (a,x: real): real; 


var 

gammcf,gln: real; 


begin 

if ( (x < 0.0) or 

begin 

writeln ( 'pause 
readln 


(a <= 0.0) ) then 
in GAMMP-invalid arguments’); 


end; 

if (x < (a+1 . 0) ) then 

begin 

gser (a,x,gammcf ,gln) ; 
gammp:= gammcf 

end 

else 

begin 

gcf (a, x, gammcf, gin) ; 
gammp : = 1.0 -gammcf 

end 

end; 


function erf {x: real) : real; 

begin 

if (x < 0 . 0) then 

begin 

erf : = -gammp (0.5, sqr (x) ) 

end 

else 

begin 

erf:= gammp (0 . 5 , sqr (x) ) 

end 

end; 
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C TWO fc pi = 6 . 283185307179586476925287 ; {value of Two pi} 
PI = 3 . 141592653589793238462643 ; {Value of PI} 
r = 6370.0; {Approximate Earth Radius} 


[GLOBAL] 

function 


INRANGE (VALUE , MIN, MAX : REAL): REAL; 


var 

NRANGES : REAL; 

N: INTEGER; 

begin __ 

if (VALUE < MIN) or (VALUE : 

begin 

VALUE := VALUE-MIN; 

MAX : = MAX -MIN; 

NRANGES := VALUE/MAX; 

N : = TRUNC (NRANGES ) ; 
if NRANGES < 0.0 then 
N : = N-l ; 

VALUE := MAX* ( NRANGES -N) ; 
INRANGE := VALUE+MIN; 

end 

else 

begin 

INRANGE := VALUE; 

end; 

end; 

[EXTERNAL] 

function MTH$ASIN(X: REAL) : 
EXTERNAL; 


- MAX) then {value outside range} 
{get how many times} 

{convert from min to max to} 

{0 to (max-min) } 

{by dividing by range} 

(get no of complete intervals} 

{make trunc work correctly} 

{get the in range value} 

{add on the minimum offset} 

{otherwise the value} 

{is already in range} 

{so do nothing} 

{back to caller} 


REAL; 


[EXTERNAL] 

function MTH$ACOS(X: REAL): REAL; 
EXTERNAL ; 


[GLOBAL] 

function 
begin 
ASIND : = 
end; 


ASIND (X : REAL): REAL; 

360 . 0 * INRANGE ( MTH$ASIN (X) , 0 , TWO_PI 


) /TWO_PI ; 


[GLOBAL] 

function SIND (X : real) : real; 

begin 

SIND : = SIN(X*TWO_PI/360 . 0) ; 

end; 

[GLOBAL] 

function COSD (X : real): real; 

begin 

COSD : = COS (X*TWO_PI/360 . 0) ; 

end; 

function angle (vll , vl2 , vl3 , v21 , v22 , v23 : real): real; 
{returns angle in degrees between two Cartesian vectors} 

var 

a,b,c: real; 

begin 

a:= vll*v21+vl2*v22+vl3*v23; 
b:= vll*vll+vl2*vl2+vl3*vl3 ; 
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c-= v21*v21+v22*v22+v23*v23 ; 

angle := MTH$ACOS (a/sqrt (b*c) ) 


*3 60. 0/TWO_PI ; 


end; 


START OP MAIN 


glndata = array [ 1 .. 1000] of real; 
qlmma = array[1..3] of real; 
glncabynca = array [ 1 . . 3 , 1 . . 3 ] of real; 
gllista = array [ 1 • • 3 ] of integer; 
glnalbynal = array [1 . .3, 1 . .3] of real, 
glcovar = array [1.. 3,1 
glnpbynp = array [ 1 . . 3 , 1 . 
glnpbymp = array [1 . .3 , 1 . 
qlnp = array [1.. 3] of integer; 
vect = array [ 1 •• 3 , 1 •• 1000] of real; 


.3] of real; 
.3] of real; 
.3] of real; 


(for MRQCOF,with 
{for COVSRT, with 
{for GAUSSJ, with 
{for GAUSSJ, with 
{for GAUSSJ, with 


nalp=3 } 
ncvm=3 } 
np=3 } 

np=3 & mp= 
np=l} 


; 3 } 


{for direction vectors) 


1 

A , DYDA : glmma; 

LISTA: gllista; 

covar, ALFA, dummy: glncabynca; 

CHISQ, ALAMDA, chisqo, deltachisq , ln0250,0500. 

Xs , XsN2 , Flux, WXs : array[1..23] of real; 

N2250 , KO , KN2 : real; . . ^ 

I, J,ndata,mf it,maxhindx,deltah: integer, 

inf i len^e^vec filename \ out f i lename : varying[80]of char; 

glochisq: real; 
glbeta: glmma; {for MRQMIN) 


{ 


-BRIGHT- 


-) 


BRIGHT 


(I: integer; H, Iinf , lnO250 ; real; var INT; real); 


procedure 

var _ _ , 

W, V, 02 50 , Vtemp, Ohm, N2hm: real ; 

hi ,h2 , zal , za2 , hm, chim : real; 
a , X , Y , F , TAU , XN2 , YN2 , HN2 , FN2 , TAUN2 ; real; 
k’, j: integer ; 

IFxWXs : real; 

posl , posm, sun : array [1.. 3] of real; 

B: real; 

begin 

0250:= exp(lnO250); 

HN2 := H*16 . 0/28 .0; 

IFxWXs := 6 . 7857e-18 ; 

Int:= 0.0; 
hi : = alt [II ; 
zal := za[I] ; 


{NB chi is SZA not SDA) 


{Integrated flux shape x xsection) 


pos 1 [ 1 ] := pv[l,I] ; 
posl [2] := pv [2,1] ; 
posl[3]:= pv[3,I]; 
sun [ 1 ] : = sv[ 1 , I] ; 
sun [ 2 ] : = sv [2,1] ; 
sun [3 ] : = sv [3 , I ] ; 
for K : = 1 to maxhindx do 
begin {integration along line of sight I) 
h2-.= hl+deltah; , . . 

za2 := ASIND ( (R+hl ) *SIND (180.0-zal)/ (R+h2) ) ; 

B : = zal-za2 ; 
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W:= SIND (B) * (R+hl ) / SIND ( za2 ) ; 
if (B = 0.0) then 
W : = h2 -hi ; 

posm [ 1 ] : = posl [l]+0.5*W*w[l,I] ; {GEI vectors for mid element} 
posm[2]:= posl [2 ] +0 . 5*W*w[2 , 1 ] ; 
posm[3]:= posl [3 ] +0 . 5*W*w[3 , 1 ] ; 

chim:= ANGLE (sun [1] ,sun[2] f sun[3] ,posm[l] ,posm[2] , posm [3] ) ; 
hm: = (hl+h2)/2; {mean altitude of element} 

if (chim >= 90.0) then {Chapman Funcs from eql7.21 of} 

{ Banks&Kockarts } 

begin 

a : = (R+hm) *C0SD {chim-90 .0) -R; {minimum ray height} 
if (a > 100.0) then 

begin 

Y : = (R+a) /H; {Y for 0} 

YN2 : = (R+a) /HN2 ; {Y for N2 } 

F:= sqrt (pi*Y/2 .0) * (1+erf (sqrt (Y/2 .0) *C0SD(chim) /SIND (chim) ) ) ; 

FN2 : =sqrt (pi*YN2/2 . 0) * (1+erf (sqrt (YN2/2 . 0) *COSD(chim) /SIND (chim) ) ) ; 
Ohm : ss 0250*exp(-(hm-250)/H) ; {[O] at hm} 

N2hm:= N2250*exp(- (hm-250) /HN2) ; {[N2] at hm} 

V : = 0.0; 

for j:= 1 to 23 do {start of wavelength loop} 

begin 

TAU : = O250*exp ( - (a-250) /H) *H*le5*Xs [ j ] *F; 

TAUN2 : = N2250*exp (- (a-250) /HN2 ) *HN2 *le5*XsN2 [ j ] *FN2 ; 

Vtemp: = Ohm*exp (-TAU) *exp (-TAUN2) *Flux[ j ] *WXs [ j ] ; 

Vt emp : = Vtemp* I inf / IFxWXs ; {normalize to parameter Iinf} 

Vtemp : = Vtemp*0. 781*0. 219/ (0 . 219+KO*Ohm+KN2*N2hm) ; 

V : = v+vtemp; 

end; {end of wavelength loop} 

INT : = INT+0 . 1*W*V ; 
end; {end of if a>100.0} 
end; {end of if chim >=90} 

if (chim < 90.0} then {Chapman Funcs from eql7.17 of } 

{Banks & Kockarts} 

begin 

X : = (R+hm) /H; 

XN2 : = (R+hm) /HN2 ; 

F : =sqrt (pi*X/2 .0) *exp ( (X/2 ) *sqr (COSD(chim) ) ) ; 

F:=F* (1-erf (sqrt (X/2 .0) *C0SD(chim) ) ) ; 

FN2 := sqrt (pi *XN2 /2 . 0 ) *exp { (XN2 /2 ) *sqr (COSD (chim) ) ) ; 

FN2 := FN2* (1-erf (sqrt (XN2/2 .0) *COSD(chim) ) ) ; 

Ohm : = 02 50* exp (- (hm-250) /H) ; { [O] at hm} 

N2hm: = N2250*exp (- (hm-250) /HN2) ; {[N21 at hm} 

V : = 0.0; 

for j:= 1 to 23 do { start of wavelength loop} 

begin 

TAU : = Ohm*H*le5*Xs [ j ] *F; {attenuation due to 0} 

TAUN2 : = N2hm*HN2 *le5*XsN2 [ j ] *FN2 ; {"* due to N2} 

Vtemp : = Ohm* exp ( -TAU) *exp (-TAUN2 ) *Flux [ j ] *WXs [ j ] ; {prod at hm} 
Vtemp := Vtemp* I inf/ IFxWXs; {normalize to parameter Iinf} 

Vtemp : = Vtemp*0 .781*0 . 2 19/ (0 . 2 19+KO*Ohm+KN2*N2hm) ; 

V:= V+Vtemp; 

end; {end of wavelength loop} 

INT : = INT+0. 1*W*V; {add contribution from chim,hm and } 

{convert to Rayleighs} 
end; {end of if chim <90} 
fit [ I ] : = INT; {keep fit to obs I for output} 
hi : = h2 ; 
zal : = za2 ; 

posl [ 1 ] : = posl [1] +W*w[l, I] ; {GEI vectors next element} 
pos 1 [ 2 ] : = posl [ 2 ] +W*w [2,1]; 
posl[3]:= posl [3 ] +W*w [ 3 , I ] ; 
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end; {end of integration along line of sight 1} 
end; {of procedure BRIGHT} 


{ 


-FUNCS- 


FUNCS 


procedure 
var 

ai,af,yi,yf: real; 

bee in 

BRIGHT { I , A [ 1 ] ,A[2] ,A[3] , y ) ; 
ai : = 0 . 999*A [ 1 ] ; 
af := 1 . 001*A[1] ; 

BRIGHT (I, ai, A[2] r A [3] ,yi) ; 
BRIGHT ( I, af ,A[2] ,A[3] ,yf) ; 
DYDA [ 1 ] : = (yi-yf ) / (ai-af ) ; 
ai : = 0 . 999 *A [ 2 ] ; 
af := 1 . 001*A[2 ] ; 

BRIGHT (I, A [ 1 ] ,ai,A[3] ,yi) ; 
BRIGHT ( I , A [ 1 ] ,af , A [ 3 ] ,yf) ; 
DYDA [ 2 ] : = (yi-yf )/ (ai-af ) ; 
ai:= 0 . 999 *A [3 ] ; 
af := 1 . 001*A[3 ] ; 

BRIGHT ( I , A[ 1 ] < A{2] ,ai,yi) ; 
BRIGHT { I , A [ 1] , A [ 2 ] ,af ,yf) ; 
DYDA [ 3 ] : = (yi-yf )/ (ai-af ) ; 
end; {end of procedure FUNCS) 


(i: integer; a: glmma; var y: real; var dyda: glmma); 


-GAUSS J- 


-} 


{ 

procedure gaussj (var a: glnpbynp; n,np: integer; var b: glnpbymp; 
integer); 

big, dum, pivmv: real; 
i , icol , irow, j , k, 1 , 11 : integer ; 
indxc, indxr, ipiv: glnp; 

begin 

for j : = 1 to n do 

begin 

ipiv [ j ] := 0 

end; 

for i:= 1 to n do 
begin 
big : = 0*0; 
for j : = 1 to n do 

begin 

if { ipiv [ j ] <> 1) then 

begin 

for k:= 1 to n do 

begin 

if { ipiv[k] = 0) then 

begin 

if (abs (a [ j , k] ) >= big) then 

begin 

big : = abs (a [ j ,k] ) ; 
irow:= j; 
icol : = k 

end 

end 

else if ( ipiv[k] > 1) then 

begin , . 

writeln ( 1 pause 1 in gauss j -singular matrix ); 

readln 

end 
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end 
end 
end ; 

ipiv[ icol] : = ipiv[icol ] +1; 
if (irow <> icol) then 

begin 

for 1 : = 1 to n do 

begin 

dum: = a [ irow, 1] ; 
a [ irow, 1] := a[icol, 1] ; 
a [icol, 1] := dum 

end; 

for 1 : = 1 to m do 

begin 

dum : = b[irow,l]; 
b[irow,l]:= b[icol,l]; 
bficol, 1] : = dum 

end 

end; 

indxr [ i] : = irow; 
indxc [i] := icol; 
if (a[icol, icol] = 0.0) then 

writeln (' pause 2 in gauss j -singular matrix ), 
readln 

end ; 

pivinv*. = 1 . 0/a t icol , icol ] ; 
a [ icol , icol ] : = 1.0; 
for 1 : = 1 to n do 

begin , . 

a [ icol , 1] : = a [ icol , 1 ] *pivinv 

end; 

for 1 : = 1 to m do 

begin . . 

b [ icol , 1 ] : = b[ icol , 1] *pivinv 

end ; 

for 11:= 1 to n do 

begin 

if (11 <> icol) then 

begin 

dum: = a [ 11 , icol] ; 
a [ 11 , icol] := 0.0; 
for 1 : = 1 to n do 

b a1n,l]:= a [11, 1] -a [icol, 1] *dum 

end; 

for 1 : = 1 to m do 

b b?n,U:= b[ll, 1] -blicol, 1] *dum 

end 

end 

end 

end; 

for 1 : = n downto 1 do 

begin 

if ( indxr [1] <> indxc [1]) then 

begin 

for k:= 1 to n do 

begin 

dum : = a [ k , indxr [111 ; 

a[k, indxr [1] ]: = a [k, indxc [ 1 ] ] ; 

a [k , indxc [ 11 ] ] : = dum 
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-} 


end 
end 
end 

end; {end of procedure gauss j } 

COVSRT 

procedure covsrt (var cover: glcovar; ncvm: integer; ma: integer; 

lista: gllista; mfit: integer), 

var 

j,i: integer; 
swap: real; 

begin 

for j:= 1 to ma-1 do 

begin 

for i:= j +1 to ma do 

begin 

covar [ i , j ] : = 0 • 0 

end 

end; 

for i: = 1 to mfit-1 do 

begin 

for j:= i+1 to mfit do 

begin 

if (lista[j] > lista[i]) then 

^ covar t lista [j 1 ,lista[i] ] := covarti, jl 

end 

else 

b covar[ lista [i] , listat jll := covarti, j] 

end 

end 

end; 

swap : = covar [1,1] ; 
for j : = 1 to ma do 

begin 

covarti, j] := covar I d,}] ; 
covar [ j , j 1 : = 0.0 

end ; 

covar [listafll ,lista[l] ] := swap; 
for j : = 2 to mfit do 

begin r - ■ 

covar [ 1 is ta[j] , lista [j] ] := covarti,]] 

end; 

for j : = 2 to ma do 

begin 

for i:= 1 to j-1 do 

begin 

covarti, jl := covar[],il 

end 

end 

end; {end of procedure covstr} 


-MRQCOF- 


•} 


{ 

procedure mrqcof tx.y.sig; glndata; ndata : integer; v.r a, glj™e ; 
ima: integer; lista; gllista; mfit; integer; ver alpha, glnalbyn , 
beta: glmma; nalp: integer; var chisq: real); 

var 

k , j , i : integer; 


var 
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ymod, wt , s ig2 i , dy : real; 

dyda : g lmma ; 

begin 

for j s= 1 to inf it do 

begin 

for k:= 1 to j do 

begin „ „ 

alpha [ j , k] : = 0.0 

end; 

beta! j 1 := 0.0 

end; 

chisq:= 0.0; 
for i : = 1 to ndata do 

begin ^ , , . 

FUNCS ( i , a , ymod , dyda ) ; 

sig2i := 1/ (siglil *sig[il ) ; 
dy : = y [ i 1 -ymod ; 
for j := 1 to mf it do 

beain . ^ , 

wt : = dyda [ lista [ j ] ] *sig2i ; 

for k : = 1 to j do 

b !lphatj,k]:= alpha! j , k] +wt*dydal lista[k] ] 

end; 

beta [ j 1 : = beta! j ] +dy*wt 

end ; . . 

chisq: = chisq+dy*dy*sig2i 

end; 

for j := 2 to mf it do 

begin 

for k : = 1 to j-l do 

begin , . . . , , 

alpha [k, j ] := alpha!;), k] 

end 
end 

end; {end of procedure MRQCOF} 

{ MRQMIN 

glncabynca; nca: integer; var chisq, alamda : real), 

label 99; 

var 

k,kk, j , ihit : integer; 
atry,da: glmma; 
oneda : glncabynca; 

begin 

if (alamda < 0.0) then 

begin 

kk : = mfit+1; 

for j : = 1 to mma do 

begin 

ihit:= 0; 

for k : = 1 to mfit do 

begin 

if (listatk] = j) then 
ihit : = ihit+1 

end; 

if (ihit = 0) then 

begin 

listatkk] := j ; 


-1 


10 



Final Report NASA Grant NAG 5-1502 APPENDIX 3 


kk: = kk+1 

end 

else if ( ihit > 1) then 

begin 

writeln ( ’pause 1 in routine MRQMIN ' ) ; 
writeln {' Improper permtuation in LI ST A' ) ; 
readln; 

end 

end; 

if (kk <> (mma+1)) then 

begin 

writeln { ’pause 2 in routine MRQMIN'); 
writeln (' Improper permtuation in LISTA'); 
readln; 

end; 

alamda:= 0.001; 

MRQCOF (x, y , sig, ndata, a, mma, lista, mf it , alpha, gibe ta, nca, chisq) ; 
glochisq:= chisq; 
for j:= 1 to mma do 

begin 

atry [ j ] := a[ j ] 

end 

end; 

for j:= 1 to mfit do 

begin 

for k : = 1 to mfit do 

begin 

covar [j ,k] : = alpha[j,k] 

end ; 

covar [ j , j ] : = alpha [ j , j ] * ( 1 . 0+alamda) ; 
oneda [ j , 1 ] : = glbeta [ j ] 

end; 

GAUSS J (covar , mfit , nca, oneda, 1,1); 
for j := 1 to mfit do 
da [ j ] : = oneda [ j , 1 ] ; 
if (alamda = 0.0) then 

begin 

COVSRT (covar , nca, mma, lista,mf it) ; 
goto 99 

end; 

for j : = 1 to mfit do 

begin 

atry [ lista [ j ] ] : = a[ lista[ j ] ] +da[ j ] 

end; 

for J : = (MFIT+1 ) to MMA do 

begin 

ATRY [LISTA [J] ] := A[LISTA[J]] 

end; 

MRQCOF (x, y , sig, ndata, atry , mma, lista, mf it , covar , da, nca, chisq) ; 
if (chisq < glochisq) then 

begin 

alamda := 0.1*alamda; 
glochisq:= chisq; 
for j : = 1 to mfit do 

begin 

for k : = 1 to mfit do 

begin 

alpha [ j , k] : = covar[j,k] 

end; 

glbeta [ j] := da[ j] ; 
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a[lista[ j] ] := atry [ lista [ j ] ] 

end 

end 

else 

begin 

alainda : = 10.0*alamda; 
chisq:= glochisq 

end; 

99: 

end; (end of procedure MRQMIN) 

GETDATA 






procedure GETDATA; 
var 

j , k , nfiles , nobs : integer; 

begin 

i : = 0 ; 

ndata : = 0 ; w - 

writeln ( 1 Enter number of input file pairs to be read 

readln (nf lies) ; 
for j:= 1 to nfiles do 

^writeln ( 'Enter name of brightness input data file: '); 

readln ( inf ilename) ; t 

writeln {' Enter name of GEI vectors input file 
readln (vecf ilename) ; 
writeln ( 'Enter # of data points : '); 

readln (nobs) ; 

open (infile, infilename, old) ; 
reset ( infile) ; 

open(vecfile,vecfilename,old) ; 

reset (vecf ile) ; 
for k:= 1 to nobs do 

begin 

ndata := ndata+1; 

readln ( infile, X [ i] ,Y[i] , sig[i] ,alt[i] , za[i] ,sza[ i] ) ; 
read {vecf ile, time [ i] , pv[ 1 , i] , pv[ 2 , i] , pv[3 , l] , 

sv 1 1 , i ] , sv [ 2 , i ] , sv [ 3 , i ] , w [ 1 , i ] , w [ 2 , 1 ] , w [ 3 , 1 ] ) ; 

ZA[ i] : = ANGLE (pv[ 1 ,i],pv[2,i],pv[3,i],w[l,i],w[2 ( i],w[3,ij) 
end; {end of one file nobs read loop} 
close (inf ile) ; 
close (vecf ile) ; 
end; {end of nfiles loop} 
writeln; , , 

writeln ( 'Enter name of fitted brightness output file • )> 

readln (outf ilename) ; 
end; {end of procedure GETDATA} 




} 


begin { main body of program} 


for I:= 1 to 3 do 
for j : = 1 to 3 do 

dummy [ i , j ]: = 0; {set dummy to zeros} 


XS[1] : = 
Xs [ 2 ] : = 
Xs [ 3 ] : = 
Xs [ 4 ] : = 


0 . 18e-18 ; {total 0 XSections 
1.3e-18; 

3 . Oe-18; 

4 . 8 e - 1 8 ; 


Richards & Torr JGR 1988} 
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XS [ 5 ] : = 
Xs [ 6 ] : = 
XS [ 7 ] : = 
Xs [ 8] : = 
Xs [9] : = 
Xs [ 10 ] : = 
XS[11] : = 
XS [ 12 ] : = 
Xs [ 13 ] : = 
Xs [ 14 ] : = 
Xs [ 15] : = 
Xs [16] : = 
Xs [ 17 ] : = 
Xs [ 18 ] : = 
Xs [ 19 ] : = 
XS [20] : = 
Xs [ 2 1 ] : = 
Xs [22 ] : = 
Xs [23 ] : = 

WXs [ 1 ] : = 
WXs [ 2 ] := 
WXs [3] : = 
WXs [4] : = 
WXs [ 5 ] : = 
WXs [6] := 
WXs [7] : = 
WXs [8] := 
WXs [9 ] := 
WXs [ 10 ] : 
WXS [11] : 
WXs [12] : 
WXs [13] : 
WXs [14] : 
WXs [15] : 
WXs [16] : 
WXs [17] : 
WXs [18] : 
WXs [19] : 
WXs [20] : 
WXs [21]: 
WXS [22] : 
WXs [23] : 

XsN2 [1] : 
XsN2 [2] : 
XsN2 [3] : 
XsN2 [4] : 
XsN2 [5] : 
XsN2 [6] : 
XsN2 [7] : 
XsN2 [8] : 
XsN2 [9] : 
XsN2 [10] 
XsN2 [11] 
XsN2 [12] 
XsN2 [ 13 ] 
XsN2 [ 14 ] 
XsN2 [15] 
XsN2 [16] 
XsN2 [17] 
XsN2 [18] 
XsN2 [19] 
XsN2 [20] 


5 . 9 e - 1 8 ; 
6 . 8e- 18 ; 
6 . 5e-18 ; 
7 . 3e-18 ; 
7 . 3e-18 ; 
8 . Oe-18 ; 
9 . le-18; 
9 . 3e-18 ; 
10 . Oe-18 
11 . Oe-18 
11 . Oe-18 
12 . Oe-18 
12 . Oe-18 
12 . Oe-18 
12 . Oe-18 
12 . Oe-18 
12 . Oe-18 
12 .Oe-18 
10 .Oe-18 


0.0373e-18;{0+(2P) R&T 
0 .276e-18; 

0 . 654e-18 ; 

1 . 08e-18 ; 

1 . 35e-18 ; 

1 . 59e-18; 

50e-18 ; 

74e-18 ; 

74e-18 ; 

01e-18 ; 

30e-18; 

2 .35e-18; 

2 . 61e-18 ; 

99e-18; 

9 6 e - 1 8 ; 

18e-18 ; 

13e-18 ; 

0 6 e - 1 8 ; 

08e-18 ; 

9 9 e - 1 8 ; 

2 ,93e-18; 

2 . 92e-18 ; 


'88 + Torr, Photchem.of Atmos. 1985} 


0 . 4 

69e- 

18; 

0 . 

60e- 

18; { 

2 . 

32e- 

18; 

5. 

40e- 

18; 

8 . 

lBe- 

18; 

9 . 

6 5e- 

18; 

10 

. 60e 

-18; 

10 

. 08e 

-18; 

11 

. 58e 

-18; 

11 

. 60e 

-18; 

14 

. 60e 

-18; 

18 

. OOe 

-18; 

17 

. 51e 

-18; 

21 

. 07e 

-18; 

21 

. 80e 

-18; 

21 

. 85e 

-18; 

24 

. 53e 

-18; 

24 

. 69e 

-18; 

23 

. 20e 

-18; 

22 

. 38e 

-18; 

23 

.10e 

-18; 


{total N2 XSections Torr Photchem. of Atmos. 1985} 
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XsN2 [21] := 23 . 20e-18 ; 

XsN2 [22] := 23 . 2 2 e - 1 8 ; 

XsN2 [23] := 29 .75e-18; 

flux [ 1 ] ; = 0.117; {Shape of modified F74113 reference spectrum} 

f lux [ 2 ] : = 0.044; 

f lux[3 ] : = 0.700; 

f lux[ 4 ] : = 0 .457 ; 

f lux[ 5] : = 0.067; 

f lux[ 6] : = 0.037; 

f lux [7 ] : = 0.257; 

f lux[ 8] : = 0.117; 

flux [ 9 ] : = 1.000; 

flux [10] := 0.143; 

fluxlll] := 0.094; 

f lux[ 12 ] := 0.047; 

f lux[ 13 ] : = 0 . 056 ; 

flux[14] := 0.041; 

flux [15] := 0.043; 

flux [16] := 0.066; 

f lux[ 17 ] := 0.104; 

f lux[ 18 ] := 0.186; 

f lux[ 19 ] := 0.051; 

flux[20] := 0.078; 

flux[21] := 0.228; 

flux[22] := 0.050; 

flux[23] := 0.033; 

GET DAT A * 

writeln(' First variable in list to be 1)H 2)Iinf or 3)[O]@250 '); 
readln(lista[l] ) ; 

writeln {' Second variable " " " 1)H 2)Iinf or 3)[O]@250 '); 

readln ( 1 ista [ 2 ] ) ; 

writeln(' Third variable " " " M 1 ) H 2)Iinf or 3)[O]@250 '}; 

readln ( 1 ista [3 ] ) ; 

writeln (' Enter the # of parameters to be adjusted: ’); 

readln (mf it) ; 

writeln ( 1 Enter initial H : '}; 

readln ( A [ 1 ] ) ; 

writeln(‘ initial Iinf : '); 

readln ( A [ 2 ] ) ; 

writelnC initial [0] @2 50 : '); 

readln ( A [3 ] ) ; 

A [ 3 ] ;= In ( A [ 3 ] ) ; 

writeln (' assumed [N2]@250: ' ); 

readln (N2250) ; 

writeln (*N2 quenching coef kN2 : '); 

readln (KN2 ) ; 

writeln (' 0 quenching coef kO : '); 

readln (KO) ; 

wr iteln (’ number of altitude intervals for BRIGHT: '); 

readln (maxhindx) ; 

writeln (' altitude interval in kilometers: '); 

readln (del tah) ; 
chisqo:= le20; 
deltachisq : = le20 ; 

alamda:= -1.0; {initialization with negative alamda) 

I : = 0; 

while (deltachisq > 0.1) do 
begin 
I : = 1+1 ; 

MRQMIN (x,y # sig, ndata, a, 3 # lista,mf it f covar # alfa, 3 , chisq, alamda) ; 
writeln ( 1 H = 1 , A [ 1 ] ) ; 

writeln{' Iinf = * , A [ 2 ] ) ; 
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writeln{' [O]@250 = ' , exp ( A [ 3 ] ) ) ; 
writeln ( ' chi sqred = ' , chisq); 
writeln ( ' end of • ,1,'th iteration 1 }; 
writeln ( 'alamda= 1 ,alamda) ; 
writeln; 

if (chisq < chisqo) then 

begin 

deltachisq : = chisqo-chisq; 
chisqo := chisq; 
end; {end of if } 
end; {end of while} 

alamda:= 0.0; {set alamda to zero for final call} 

MRQMIN (x,y , sig, ndata, a, 3 , lista, mf it , covar , alfa, 3 , chisq, alamda) ; 

writeln ( 'alpha matrix:-'); 

writeln (alfa [ 1 , 1] , alfa [1,2] , alfa [1,3] ) ; 

writeln (alfa [ 2 , 1 ] , alfa [ 2 , 2 ] , al fa [2 , 3 ] ) ; 

writeln (alfa [3 , 1] , alfa [3 , 2 ] , al fa [3 , 3 ] ) ; 

writeln; 

writeln (' covariance matrix:-'); 
writeln (covar [1,1] , covar [1,2] , covar [1,3]); 
writeln (covar [2,1], covar [2,2], covar [2,3]); 
writeln (covar [3,1] , covar [3 , 2 ] , covar [3 , 3 ) ) ; 
writeln; 

writeln ( * H = ' , A [ 1] s 7 , ' +/- 1 , sqrt (covar [ 1 , 1] ) : 7); 

writeln { ' I* = 1 , A [ 2 ] : 7 , ' +/- 1 , qrt (covar [2 , 2 ] ) : 7); 

writeln ( ' [0] @250 = ',exp(A[3]): 7 100* (exp (sqrt (covar [3 , 3 ])) - 
1.0) : 7, ■%'); 

0500 : = exp (A [3 ] ) *exp ( (250. 0-500. 0 ) / A [ 1 ] ) ; 
writeln{' [0]@500 = ',0500: 7); 

0500 : = exp (A [3 ] ) *exp ((250.0-500.0)/ ( A [ 1 ] +sqrt (covar [1, 1] ) ) ) ; 
writeln{' [0]@500+= ',0500: 7); 

0500 := exp (A [3 ] ) *exp ((250.0-500.0)/ (A [ 1 ] -sqrt (covar [1 , 1] ) ) ) ; 
writelnC [0]0500-= ',0500: 7) ;writeln; 
writeln ('chi sqred = ',chisq);writeln; 
writeln; 

open (outfile, outfilename, new) ; 
rewrite (out file) ; 
for i : = 1 to ndata do 

begin 

write (outfile, X[i] , ' ',fit[i],' ',y[i], r *); 
writeln(outfile,sig[i],' ' , za[ i ] , ' SZA[ i ] ) ; 
end; {end of read loop} 
close (outfile) ; 

end. {end of program TWIFITTER} 
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